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Abstract. The aim of this review is to provide an overview of a recent work 
concerning "leaky" quantum graphs described by Hamiltonians given formally 
by the expression —A — a5{x — F) with a singular attractive interaction sup- 
ported by a graph-like set in R", u = 2,3. We will explain how such singu- 
lar Schrodinger operators can be properly defined for different codimensions 
of r. Furthermore, we are going to discuss their properties, in particular, 
the way in which the geometry of F influences their spectra and the scatter- 
ing, strong-coupling asymptotic behavior, and a discrete counterpart to leaky- 
graph Hamiltonians using point interactions. The subject cannot be regarded 
as closed at present, and we will add a list of open problems hoping that the 
reader will take some of them as a challenge. 
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1. Introduction 

In this paper we are going to review results concerning a class of "different" 
quantum graph models. With this aim in mind, it would be natural to start by 
recalling briefly the standard quantum graphs, their description, properties, and 
numerous applications. In this volume, however, this would be clearly to bring 
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owls to Athen|3 and we refrain from doing that referring to the other articles in 
these proceedings, or to |BCFK06] as another rich bibliography source. 

To motivate a need to look for alternative description of graph-like structures, 
let us observe that — despite its mathematical simplicity, beauty, and versatility 
— the standard quantum-graph model has also some drawbacks. In our opinion, 
the following two are the most important: 

• the presence of ad hoc parameters in the boundary conditions which de- 
scribe how the wave functions are matched at the graphs vertices 

• the fact that particles are strictly confined to graph edges. While this 
is often a reasonable approximation when dealing, say, with electrons in 
semiconductor quantum wires, such a model neglects quantum tunneling 
which can play role once such wires are placed close to each other. One 
consequence is that in such a description, with the graph embedded in W^, 
spectral properties reflect the topology while the graph geometry enters 
only through the edge lengths, their shapes being irrelevant 

A way to deal with the first problem is to regard a quantum graph as an idealization 
of a more "realistic" system without such ambiguities; a natural candidate for this 
role are various "fat graphs" . Limits of such objects when the edge width squeezes 
to zero were studied extensively, first in the easier Neumann-type case |FW93|, 
KZOl, RSOl, SaOl, Ku02, KZ03, EPOS, Po06, EPOr and quite recently also 
in the situation with Dirichlet boundaries [PoOSl IMV071 ICE071 IGrOT) . These 
results give a partial answer to the first quest ior0 while the second problem remains. 

Here we are going to discuss a class of quantum graph models which are free 
of both difficulties; the price we pay is — similarly as for the fat graphs mentioned 
above — that instead of ordinary differential equations we have to deal with a PDE 
problem. The idea is to preserve the whole Euclidean space as the configuration 
space and to suppose that the particle is kept in the vicinity of the graph F C K."^ 
by an attractive singular interaction. Formally such a Hamiltonian expresses as 

(1.1) - A - aix)6{x -T) 

with a{x) > ; we will consider mostly the situation where the attraction is position 
independent, a{x) — a > 0. Before proceeding to a definition of such singular 
Schrodinger operators and discussion of their properties, let us make a few remarks. 

First of all, it is clear that there is no ambiguity related to the graph vertices 
once F and a are given. It is equally obvious that the confinement in this model 
takes place at negative energies only. The particle now "lives" in the whole space 
and can be found even at large distances from F, although with a small probability, 
because the complement M'' \ F is the classically forbidden region. The presence of 
the tunneling is the reason why we dub such systems as leaky quantum graphs. 

Schrodinger operators with interactions supported by curves and other mani- 
folds of a lower dimension were studied already in the early nineties |BT92j . and 
even earlier in examples with a particular symmetry [AGSW, ISha88| . A more 
systematic investigation motivated by the above considerations was undertaken in 
a last few years; it is the aim of this review to describe its results. 



^Or, depending on your taste, coal to Newcastle, firewood to the forest, etc. As usual, one 
can also refer to the Bard: to throw a perfume on the violet (The Life and Death of King John). 

^An alternative approach is to keep the graph fixed and to approximate the vertex coupling 
through suitably scaled families of regular or singular interactions - see [Ex96L ICE04L lETuOT] . 
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One should stress, however, that such mathematical structures can be studied 
also from other points of view. A prominent example comes from studies of high 
contrast optical systemsQ used to model photonic crystals — see, e.g. |FKu96l 
IKK98| — which in a suitable approximation yield an analogue of the spectral 
problem for the operator (jl.ip : the two differ only by the physical interpretation, 
the roles of the coupling and spectral parameters being switched. A derivation of 
leaky-graph models in this context was given in the paper FKu98], see also the 
review |Ku01| and recall that the corresponding operators can be cast also in a 
pseudo-differential form TKu98l [PP04] . 

The material we are going to review is relatively extensive. We will take care, 
of course, to explain properly all the notions and the results. On the other hand, 
proofs will be mostly sketched. However, we will always give references to original 
papers where the particular complete argument can be found. Let us finally remark 
that the subject reviewed here cannot be regarded as closed, on the contrary, there 
are many open questions. We devote to them the closing section, and the author of 
this survey can only hope that his reader will take this problem list as a challenge 
and a program which will keep him or her busy for some time. 



2. Leaky graph Hamiltonians 

2.1. Quadratic forms and boundary conditions. The Hamiltonians we 
are interested in arc generalized Schrodingcr operators with a singular interaction 
supported by a graph- like F which is a zero me asure set i n M.'^. We will use facts 



about such operators derived, in particular, in |BEKS94j specifying them to our 
present purpose. Let us first suppose that the configuration space dimension 1^ — 2 
and the coupling "strength" is constant on the interaction support. 

To begin with, let us show how such a singular operator can be defined generally 
through the associated quadratic form. Consider a positive Radon measure m on 
and a number a > such that 



(2.1) (1 + a) / \Tp{x)\^ dm{x) < a \\'ip{x)\^ dx + b \ip{x)\^ dx 

JR2 JR2 JM2 

holds for all ip € 5(E^) and some a < 1 and b. The map defined by I„iip = ?/; 
on the Schwartz space 5(R^) extends by density uniquely to 

(2.2) : M^1'2(E2) ^ L2(^rri) := L^{R'^,m) ; 

for brevity the same symbol is used for a continuous function and the corresponding 
equivalence classes in both L^(M'^) and L^(m). The inequality (|2.ip extends to 
^i,2^]]j2-j ^ replaced by the left-hand side. The quadratic form 

(2.3) S-a^mm-.^ [ \\7iP{x)\''dx-a [ \ilm^){x)\^ dm{x) 



is defined on W ' (R ); it is straightforward to check [BEKS94| that under the 



condition (|2.1[) this form is closed and below bounded, with C^(R ) as a core, and 
consequently, it is associated with a unique self-adjoint operator denoted as H-am- 



Another situation where one arrives at a leaky-graph-type model arises when one deals with 
contact interactions of several one-dimensional particles jDuOS,. 
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A sufficient condition for tlic inequality (12. 1[) to be valid is that the measure m 
belongs to the generahzed Kato class, i.e. 

(2.4) lini sup / I In jx — ?/|| dm(y) = , 

where B^(x) is the ball of radius e centered at x. In such a case, moreover, any 
positive number can be chosen as a. So far the construction has been general 
and involved also regular Schrodinger operators. Suppose now that m is the Dirac 
measure supported by a graph F C which has the following properties: 

(gl) edge smoothness: each edge ej € F is a graph of function jj : 
Ij M.'^ where Ij is an interval (finite, semi-infinite, or the whole R). 
Moreover, without loss of generality we may suppose that edges are para- 
metrized by the arc length, |7j(s)| = 1. 

(g2) cusp absence: at the vertices of F the edges meet at nonzero angles. 

(gS) local finiteness: each compact subset of contains at most a finite 
number of edges and vertices of F. 



The last assumption allows us to extend Theorem 4.1 of [BEKS94| applying it to 



the Dirac measure supported by the graph. More exactly, we consider the measure 

(2.5) mr : mr(Af) = £i{M n F) 

for any Borel Af C K'^, where £i is the one-dimensional Hausdorff measure given 
in our case by the edge-arc length. Such a straightforward extension implies that 
m = mr satisfies the condition (j2.ip and gives thus rise to the appropriate operator 
H^am] to make the notation explicit we will employ for it in the following the 
symbol Hax- This is one way how to give meaning to the formal expression (jl.ip . 
An alternative is to use boundary conditions. Consider the operator acting as 

(7i„,rV) (a;) --(AV)(a;), x G \ F , 

on any function ip which belongs to VK^'^(R^ \ T), is continuous at each edge ej E F 
with the normal derivatives having there a jump, namely 

dtp dtp 

(2.6) -^^{x) - -^^{x) = -atp{x) , xe'mtej; 

since the edges are smooth by assumption, the norma l vector exists at each inner 
point of an edge. In the same way as in [BEKSoS] one can check that Ha,r is 
e.s.a., and moreover, by Green's formula it reproduces the form ()2.3p on its core, 
so its closure may be identified with Hax defined above. 

Remarks 2.1. (i) The above definitions easily extend to the situation with the 
singular interaction strength a{s) varying along the edges provided the correspond- 
ing function a : F — > M_|_ is sufficiently regular. 

(ii) In a similar way one can define operators corresponding to the formal expression 
(jl.ip for a generalized "graph" whose edges are {ly — 1) -dimensional manifolds in 
R'' satisfying suitable regularity conditions analogous to (gl)-(g3). 

2.2. Regular potential approximation. As we shall see below the opera- 
tors Ha^r represent a reasonably general class of systems for which various prop- 
erties can be derived. One can ask nevertheless whether this is not again a too 
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idealized model. Before proceeding further we want to show that Ha^r can be re- 
garded as weak-coupling approximation to a class of regular Schrodinger operators; 
for simplicity we restrict ourself to graphs with a single infinite edge. 

Let r be a curve described by a function 7 : M ^ M^. Then we are 
able to define the signed curvature k(s) := (7172 — 7271) (s); we assume that it is 
bounded, |fc(s)| < c+ for some > and all s e M. Moreover, we suppose that 
r has neither self-intersections nor "near-intersections" , i.e. that there is a c_ > 
such that |7(s) — 7(s')l ^ f'^'' ^^V/ 5,3' with |s — s'| > c_. Then we ca n defin e 
in the vicinity of F the standard locally orthogonal system of coordinates ES89J , 
i.e. the pairs (s, u) where u is the (signed) normal distance from T and s is the 
arc-length coordinate of the point of F where the normal n(s) is taken; the system 
is unique in the strip neighborhood :— {x{s, u) : (s, u) g E^}, where 

(2.7) a;(s, u) := 7(5) + n(s)u 

and E° := {{s,u) : s G M, |u| < e} as long as the condition 2e < c_ is valid. 

With these prerequisites we can construct the approximating family. Given 
W G L°°((— 1, 1)), we define for all e < ^ c_ the transversally scaled potential. 



(2.8) V,{x) 



... X^Ee 

... ^eE, 

and put 

(2.9) H,{W,-f):^-A + V,. 

The operators iJe(VF, 7) are obviously self-adjoint on D(— A) = T/F^'^(R^) and we 
have the following approximation result: 

Theorem 2.2. Under the stated assumptions, H^{W,T) — > H^.r as e — > in 
the norm-resolvent sense, where a J^-^^ W{t)dt. 

Sketch of the proof. One has to compare the resolvents, that of 
given below and the Birman-Schwinger expression of {H,{W,-f)-k^) . Both are 
explicit integral operators and their difference can be treated in a way similar to 
that used in the squeezing approximation of the one-dimensional S interaction — 
see, e.g., [AGHH04) — a full account of the argument can be found in |EI01] . □ 

Notice that the regular potential approximation is not the only way how to 
justify the leaky-graph model physically; for an alternative see [FKu98j. 

2.3. The resolvent. As usual, the spectral and scattering properties are en- 
coded in the resolvent and our first task is to find an explicit expression for this 
operator. We will employ an analogue of the Birman-Schwinger formula for our sin- 
gular case. If fc^ belongs to the resolvent set of i?a,r we put i?^ p ■— {Hax — k^)~^- 
We look for the difference of this operator and the free resolvent Rq which is for 
Im fc > an integral operator with the kernel 



(2.10) Gk{x^y) = ^Hl,'\k\x-y\). 



4 



To this aim we need embedding operators associated with Rq. Let fi, v be arbitrary 
positive Radon measures on with ^{x) = v{x) =0 for any x £ M^. By i?^ we 
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denote the integral operator from L^(^) :— i^(M^, fi) to L?{v) with the kernel Gk, 
i.e. 

holds i^-a.e. for all (f) E D{R'l ^) C L'^ilj)- In our case the two measures will 
be m = mr introduced by (|2.5p and the Lebesgue measure dx on in different 
combinations, which simply means that one or both variables in the kernel ()2.10p 
are restricted to T. Using this notation we can state the following result: 

Proposition 2.3. (i) There is a kq > such that the operator I — ald-^ ^ on 
iP'im) has a bounded inverse for any k> kq. 

(ii) Let Im/c > 0. Suppose that I — aR^-^ ,„ is invertible and the operator 
R'^ := Rq + ai?j2.^[/ — Q;i?f„ ,„] ^Rm,dx 

from L^(]R.^) to L^(M^) is everywhere defined. Then belongs to p{Har) and 
(i/„,r-P)-i = i?'=. 

(Hi) dimker(i?Q_r — = dimker(/ — aiij^j for any k with Imfc > 0. 

(iv) an eigenfunction of Ha.v associated with such an eigenvalue fc^ expresses as 

ip{x) = i?^^_,^(x,s)(/)(s)ds, 
Jo 

where cj) is the corresponding eigenfunction of aR^ „ with the eigenvalue one. 

Sketch of the proof. The result, which is in fact valid for any operator 
H-am of Sec. 12. H is obtained by verifying the Birman-Schwinger (BS) formula in 
our singular setting. The procedure requires some care; a full account concerning 
the claims (i)-(iii) can be found in BEKS94 , for (iv) see |Pos04j . □ 

2.4. The case of codimension two. The second one of the above construc- 
tions of Ha.r can also be rephrased in the following way: we first restrict the 
Laplacian to a symmetric operator defined on function which vanish in the vicinity 
of r, and afterwards we choose a particular self-adjoint extension specified by the 
condition l|2.6p . It follows from general properties of partial differential operators 
|He89j that a similar construction is possible also in higher dimensions as long as 
codimF < 3. Since we want to stick to cases of physical interest, we will mention 
here only graphs whose edges are curves in R^. 

An analogue of the form definition (|2.3p does not work in this situation and 
we have to rely on boundary conditions which are, however, more complicated 
than (|2.6p . The difference is of a local character, thus we restrict ourself to the 
simplest situation when F is a single infinite curve described by a function 
7(s) : R ^ without self-intersections and such that |7(s)| = 1. In view of the 
smoothness assumption the curve possesses locally Frenet's frame, i.e. the triple 
(t(s), b{s),n{s)) of the tangent, binormal and normal vectors; we assume its global 
existencqj. For a fixed nonzero p G R.^ we define the "shifted" curve Fp as the 
graph of the function 

7p(s)) := 7(s) + Pib{s) + P2n{s) ; 



This is true, for instance, if 7 does not vanish. For curves having isolated straight segments 
a suitable coordinate system can be obtained by patching local Frenet systems together, possibly 
with a rotation - cf. a discussion on that point in |EK04a . 
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the distance between the two is obviously r Ipj. If we suppose in addition that 
r does not have "near-intersections" as in Sec. 12.21 then clearly Fp n F = holds 
provided r is small enough. Since any function / S Wf^'^{M? \ F) is continuous on 
R'^ \ F its restriction to Fp is then well defined; we denote it as /L , in fact, we 

can regard /L as a distribution from D'{R) with the parameter p. We denote by 

V the set of functions / e Wf;l{E!^ \ T) n L'^{M!^) such that the following limits 



r^o m r p 
n{f){s) := hm/L (s) + S(/)(s)lnr 



+0 



exist a.e. in R, are independent of the direction ip, and define functions from 
L^(R); the limits here are understood in the sense of the Z3'(R) topology. Now we 
are able to define the singular Schrodinger operator in the present case: it acts as 

(2.11) iJ„,r/ = -A/ for a; e R^ \ F 
on the domain 

(2.12) D{H^s) ■■={geV: 27:aE{g){s) = n{g)is)} . 

In this way we get a well-defined Hamiltonian which we seek |EK02| : 

Theorem 2.4. Under the stated assumptions the operator Haj- is self-adjoint. 



As in the case of codimension one and the conditions (|2.6p the above definition 
has a natural meaning as a point interaction in the normal plane to F. 

The proof of Theorem l2.4l is technically more involved and we restrict ourselves 
to a few remarks referring to |EK02j for the full exposition. The argument is in a 
sense opposite to the previous considerations. It is based on an abstract analogue 
of Proposition 12.31 proved in |Pos01| , see also |Pos04| , which shows existence of 
a self-adjoint operator with the resolvent of the appropriate form, after that one 
verifies that this operator coincides with the above Ha.r- 

Remark 2.5. Since the mentioned resolvent formula analogous to that of 
Proposition 12.31 will be used in the following, we will describe it at least briefly. 
It contains again traces of the free resolvent, which is now given by 

(2.13) Gkix-y)^— 

An \x — y\ 

However, the use of Posilicano's abstract result requires to interpret the embedding 
operators involved not as maps between spaces, but rather the last factor i?p as 
M^2,2(]^3) _^ 2.2(R) and its counterpart as the Banach space dual to i?p. Then 

(2.14) R'' ^ RI^~ R'^{Q''--a)-^R'^, 

where the modified position of a in this formula corresponds to the usual convention 
about the coupling parameter for two-dimensional point interactions [AGH H04] 
reflected in (|2.12p : roughly speaking, it is an inverse of the one appearing in p.6p . 
The operator Q'' is the counterpart to „ of Proposition 12.31 but we use on 
purpose a different symbol to stress that we cannot write it simply as an integral 
operator and a renormalization is needed, cf. |EK02j for more details. 
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3. Geometrically induced properties 

We said in the opening that even if we think about F in the usual quantum 
graph model as embedded in R'', the shapes of the edges do not influence the 
spectrum. Leaky graphs are different as one can illustrate in various ways. 

3.1. Bound states due to non-straightness. Consider again a leaky graph 
in M^. If F = Fq is a straight line corresponding to 70 (s) = as + b for some a,b ^M.^ 
with \a\ = 1, we can separate variables and show that 

1 



■ , 00 



(3.1) ^(i?a,.„. L 4 

is purely absolutely continuous. We are going to show that a bend or deformation 
produces, within a wide class of curves F, a non-void discrete spectrum. To be 
specific, we assume that the generating function 7 : R — > is continuous and 
piecewise (or, in terms of the assumption (gl)-(g3), a graph which may have 
vertices but no branchings) satisfying the following condition^j: 

(al) a lower distance bound: there is c € (0, 1) such that |7(s) — 7(s')| > 

c\s — s'\. In particular, F has no cusps and self-intersections, and its 

possible asymptotes are not parallel to each other. 
(a2) asymptotic straightness: there are positive d, fi > ^, and uj E (0, 1) 

such that the inequality 

^_ hlf)^2p<d[i + |, + ,'|2.]-V2 

holds true in the sector S^) ■= {{s, s') : oj < < uj^^ }. 
(a3) non-triviality: we excluded the case F = Fq. Recall that 

|7(s) -7(s')| < \s- s'\ 

holds for any s, s' E M, hence we request in other words that the last 
inequality is sharp at least for some s, s' S M. 
Then we have the following result: 

Theorem 3.1. Let a > and suppose that 7 : R ^ satisfies the above 
assumptions. Then the essential spectrum is the same as for the straight line, 
o'css{Ha.r) — but Ha.r has at least one isolated eigenvalue below —jo^. 

Sketch of the proof. Observe first that in view of (a2) it is not difficult to 
construct a Weyl sequence to Ha^r showing that any non-negative number belongs 
to (Tess- To deal with the negative part, we use the generalized BS principle of 
Proposition 12.31 the idea is to treat the difference between the operator TZ'^ p := 
m on i^(R) and its counterpart corresponding to Fq as a perturbation. The 
integral kernel of the operator is 

7^S,^(5,s')-£^oW7W-7(s')l) , 

where ifo is the Macdonald function; for F = Fq one has to replace |7(s) — 7(s')l 
by I s — s' I . In the last named case the operator is of convolution type and using 
Fourier transformation it is easy to check that its spectrum is absolutely continuous 
covering the interval [0, Q!/2k] in correspondence with ()3.1|) . 



If 7 S we have a sufficient condition for (a2) in terms of the signed curvature introduced 
Sec. [22] it is valid with > 5 if k{s) = 0(\s\-f) with /3 > | as |s| ^ 00, of. | EI01| . 
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The key observation is that the kernel of :— TZ^^r ~ ^q.Fo satisfies 

(3.2) v,{s,s') := ^ (^Ko{k\^{s)-j{s')\) ~ Ko{n\s-s'\)^ >0 

in view of (al) and the monotonicity of Kq, the inequahty being strict for at least 
some values of the variables s, s'. We shall then argue in three steps: 

(1) a variational argument in combination with (|3.2p shows that the spectrum 
is "pushed up" by the perturbation, sup cr (jZ'^ p) > ^ if F is not straight. 

(2) in view of (a2), 2?^ is Hilbert-Schmidt for fi > ^ , and therefore compact. 

(3) the map k i-^ TZ'^ p is operator-norm continuous and TZ^.-y ^ as k ^ oo. 
The compactness of implies, in particular, in combination with Proposition l2.3l 
the claim about the negative part of the essential spectrum. 

The discrete spectrum part can be also derived from the generalized BS prin- 
ciple. It follows from the above claims that there are spectral points of TZ^. p above 
a /2k and they cannot be anything but eigenvalues of a finite multiplicity. More- 
over, every such eigenvalue depends continuously on k and tends to zero as k ^ oo. 
Hence it crosses one at a value kq > giving rise to the sought eigenvalue of the 
operator Ha^v', for details of the argument see [EIOl] . □ 

It may seem that the result covers only a rather particular class of graphs. 
Using the minimax principle, however, we arrive at the following easy consequence: 

Corollary 3.2. Suppose that T has a subgraph in the form of an infinite curve 
satisfying the assumptions of the theorem, and acss{Ha,T) — then the 

discrete spectrum ofHa Y is non-empty. 

It is important that the assumption about preservation of the essential spectrum 
can be often verified easily, for instance, in the situation when F has outside a 
compact a finite number of straight edges separated by non-trivial wedges. 

3.2. An example: leaky star graphs. To illustrate the last made claim, let 
us investigate in more detail a particular class of such graphs, namely the situation 
when F is of a star shape. Given an integer N > 2^ consider an (A'^— l)-tuple 
[3 = . . . , /Stv-i} of positive numbers such that 

Pn := 27r - ^ /3, > . 

j=i 

Denote :— ^^"^ '■= 0. Let Lj be the radial half-line with the endpoint 

at the origin, Lj := {x G ]S? : argx = -dj}, naturally parametrized by its arc 
length s — \x\. These half-lines will be the edges of F = F^ := UjlJo^ to which 
we associat^ the Hamiltonian Hn{(3) ■— Ha^Vf,- Trivial examples are 

(1) H2(tt) corresponding to a straight line which has obviously a purely a.c. 
spectrum covering the interval [— ia^,oo), 

(2) H4{Ps) with f3s = {f : f > f } corresponding to cross-shaped F allows again 
a separation of variables. The a.c. part of its spectrum is the same as 
above, and in addition, there is a single isolated eigenvalue —^o? corre- 
sponding to the eigcnfunction (2a)~^e~"^l^l"'"l^l^/^. 



Properties of Hiq{f3) depend on the order of the angles, however, operators related by a 
cyclic permutation are unitarily equivalent by an appropriate rotation of the plane. 
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Star-shaped graphs have the property indicated above: 

Proposition 3.3. acasiHpf{P)) — [—ja'^,00) holds for any N and (3. 

Sketch of the proof. The inequahty inf tTcss(-ffAr(/3)) > is obtained 

by Neumann bracketing dissecting the plane outside a compact into semi-infinite 
strips with Lj in the middle and "empty" wedges. The fact that [— |;a^, 00) belongs 
to the spectrum is checked by means of a family of Weyl sequences, cf. |EN03| . □ 

By Corollarv l3.2i (Tdisc(^Jv(/3)) is nonempty unless N = 2 and /3 = tt. Using 
direct methods one can establish various other properties of the discrete spectruixlj. 

Proposition 3.4. Fix N and a positive integer n. If at least one of the angles 
[3j is small enough, ft crdisc(-ffAf(/3)) > n. In particular, the number of bound states 
can exceed any fixed integer for N large enough. 

Sketch of the proof. By the minimax principle it is sufficient to check the 
claim for H2{l3). We choose the coordinate system in such a way that the two 
"arms" correspond to aigO — ±/3/2 and employ trial functions of the form y) — 
f{x)g{y) supported in the strip L < x < 2L, with f £ satisfying f{L) = f{2L) = 
0, and 

1 ... |y|<2d 

^-a{\y\-2d) \y\>2d 

with d := Ltan(/3/2). Evaluating the quadratic form of i?2(/3) and using minimax 
principle again we get the result, see [EN03 for details. □ 



giy) = 



Remark 3.5. The above variational estimate also shows that the bound state 
number for a sharply broken line is roughly proportional to the inverse angle, 

o3/2 

as /? — *■ 0. This can be regarded as an expected result, since the number is given by 
the length of the effective potential well which exists in the region where the two 
lines are so close that they roughly double the depth of the transverse well. 

Let us see how the BS equation looks like explicitly for star graphs. Define 



(3.3) d^j{s,s') = d^j{s,s') = •y/s2+s'2-2ss'cos|t?j - 



with - = EL+i A, in particular, du{s,s') = \s - s'\. By 7^« (/3) = 7^^,(/3) 
we denote the operator L^(K+) L2(M+) with the integral kernel 

nfj{s,s';P) — i^o(Kd^,(s,s')) ; 

then the (discrete part of the) spectral problem for the operator IIn{(3) is by Propo- 
sition 12.31 equivalent to the matrix integral-operator equation 

N 

(3.4) ^(7^r^.(/3)-5.,/)0, =0, z = l,...,iV, 

j=i 



^Variational methods can be also used to establish the existence of discrete spectrum in some 
range of parameters, see the paper [BEW 08 in this volume. 
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on L^(M+). Notice that the "entries" of the above kernel have a monotonicity 

property, 7^« (/?) > TZ^jifi') if - < I'&'j - This fact has the foUowing easy 
consequence |EN03j : 

Proposition 3.6. Each isolated eigenvalue A„(/3) of H2{(3) is an increasing 
function of the angle [3 between the two half-lines in (0,7r). 

3.3. Higher dimensions. If we restrict ourselves to the physically interesting 
case of three-dimensional configuration space, there are two possible ways how 
to extend the above results. One refers to the situation when the interaction is 
supported by a surface. Here unfortunately only a particular result is known to the 
date, which we will mention in Sec. 14.21 below. 

Let us thus consider the second possibility when F is an infinite piecewise 
curve in M^. The argument is similar to that of the previous section but it 
needs care due to the more singular character of the interaction. If F = Fq the 
spectrum of i?Q,r is found by separation of variables; the known properties of the 
two-dimensional point interaction [AGHH04j imply that 

where Cq = -4e'^(-'^-^"+^W^- in this expression « 0.577 is Euler-Mascheroni 

constant. To rephrase it in terms of the BS operator Qq :— Qq^ notice that the 
latter equals 7^ + ln2 + 7/'(l) where 7^ in the momentum representation acts as 
multiphcatior0 by ln(p^ -I- k^)^/^. Consequently, the spectrum of Qq is purely 
absolutely continuous and equal to (—00, s^] where := ^ii^i^) ~ lii('*/2))- 

To proceed we need assumptions about the curve F. We retain (al) and (a3) 
from the previous section, while (a2) will be replaced by 

(a2') ASLS: there are positive d, fi > ^, and uj E (0, 1) such that 
. l7(^)- 7(^01 ^ , k^£l 

|S-S'| - (1+ |s-s'|)(l + (sHs'2)A')l/2 

holds true in the sector Si^, the same as before. 

In addition to the asymptotic straightnes^ we require newly some local smoothness 
of the curve. Now we can follow proof of Theorem 13. II step by step; after checking 
that the essential spectrum is preserved, and denoting Q'^ :— Q***, we prove that 

(1) sup cr(Q'^) > Sk by a variational argument using the sign definitess, 

V,{s,s') ^ G,«(7(s)-7(s')) - G,«(s-s') > 0, 

with a sharp inequlity at least for some values of the variables s, s'. 

(2) in view of (a2') is Hilbert-Schmidt for /i > ^ , hence compact, and 

1/2 

the corresponding norm ||2?K|jjjg is uniformly bounded w.r.t. k > \(a\ 

1/2 

(3) the function k is operator-norm continuous in (|Ca| ,00) and 
Q'^ —00 holds as K ^ 00. 

Working out details of this scheme [EK02) we arrive at the following conclusion: 



'^Comparing to {27r)^''^ (p'^ + k?)^^^"^ in the codimension one case. This is why the kernel of 
Qg makes sense as a distribution only and a renormalization mentioned in Remark l2.5l is needed. 

^For C'^ smooth curves we have a sufficient condition for the appropriate large-distance 
behavior analogous to that mentioned in the footnote to assumption (a2). 
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Theorem 3.7. Fix a > and suppose that the generating function 7 ; R — > R'^ 
0/ r satisfies the stated assumptions. Then CToss(^fa,r) ~ [Ca,oo), and the operator 
Ha,T has at least one isolated eigenvalue in the interval (— cx),<^q). 

Let us remark that the strengthened hypothesis in (a2') is not needed to prove 
the existence of the geometrically induced spectrum, but rather to determine its 
character; due to the more strongly singular character of the interaction in the 
codimension two case the local smoothness is required to guarantee its discreteness. 

3.4. Geometric perturbations. Let us turn to another way in which the 
leaky character of our graphs is manifested. Consider a graph F C M'' with two 
edges, the endpoints of which are close to each other; we can think of this situation 
as of a single edge having a hiatus. In the standard quantum graph setting, it is 
only the topology which matters, either the two edges are connected or not. Here, 
in contrast, the distance of the two endpoints plays a role. 

Let us first analyze the general codimension one situation, the (/-dimensional 
Schrodinger operators with a ^-interaction supported by a (1/ — l)-dimensional 
smooth surface having a "puncture" . We are going to show, formally speaking, that 
up to an error term the eigenvalue shift resulting from removing an e-neighborhood 
of a surface point is the same as that of adding a repulsive 5 interaction at this 
point, with the coupling constant proportional to the puncture "area". 

Let F be a compact C"'-smooth surface in W with r > ^v; without loss of 
generality we may suppose that it contains the origin. Let further {7'e}e>o be a 
family of subsets of F which obeys the following requirements: 

(pi ) measurability: Ve is measurable with respect to the (i/— l)-dimensio- 

nal Lebesgue measure on F for any e small enough. 
(p2) shrinking: sup^g^^ \x\ = 0{e) as e — > 0. 

Consider the operators Ha,r and Ha^r, corresponding to F^ :—T\V(, defined as in 



Sec. 12.11 Since F^ is bounded, we have 

acssiHaXe) = [0,00) and tt o-disc(i?Q.r) < 00 . 

By the minimax principle there is a unique a* > such that (Tdisc(^ct,r) is non- 
empty if a > a* while the reverse is true for a < a*; it is not difficult to check that 
when = 2 and a* > for j/ > 3, see jBEKS94l Thm 4.2]. 



-■ 

Let TV be the number of negative eigenvalues of Ha,T- Using the convergence of 
the corresponding quadratic forms (|2.3p on W^''^{E}') as e ^ we find by [Ka76| 
Thm VIIL3.15] that for all e small enough Hax, has the same number N of negative 
eigenvalues, which we denote as Ai(e) < A2(e) < • • • < AAr(e), and moreover 

Aj(e) ^ Aj(0) as e -> for l<j<N. 

Let {'■Pj{x)}^^i be an orthonormal system of eigenfunctions of H^.t corresponding 
to these eigenvalues; without loss of generality we may suppose that ipi [x) > 
in R"^. Using the Sobolev trace theorem, one can check that each function ipj is 
continuous on a F-neighborhood of the origin. Given S UdisciHa^) we denote 
as and n{i') the smallest and largest index value j, respectively, for which 

/i = Aj(0), and we introduce the positive matrix 



(3.5) 
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denoting by Sm(/i) ^ Sm{ii)+i < • • • < ^n{ti) its eigenvalues. In particular, if /i = 
Aj(0) is a simple eigenvalue of Ha.r, we have m{fi) = n{y) = j and Sj = |(/7j(0)p. 
With these prerequisites, we can make the following claim: 

Theorem 3.8. Assume (pi), (p2), and suppose that a > a* . For a given 
^ (z o'disciHa,r) we have the asymptotic formula 

Xj{e) — fj, + amY{Ve)sj + o{e'^^^) , m{fi) < j < n{^) , as e^O, 

where r7ir(') stands for the {v— I) -dimensional Lebesgue measure on T. 

To derive the expansion, one has to note that due to the singular character of 
the perturbation a direct use of the asymptotic perturbation theory of quadratic 
forms |Ka76[ Sec. VIII. 4] is not possible. Indeed, we have 

for V e C^iR") and the quadratic form Cg°{W) 9 V ^ IVXO)^ £ R does not 
extend to a bounded form on W^''^{M.'^), because the set oi ip E C^(R'^) vanishing 
at the origin is dense in M^^'^(K'^). A way to eliminate this difficulty is to employ 
the compactness of the map W'^''^{R'')3 f ^ f\r ^ ^^(r); we refer to [EY03j for 
a detailed description of such a proof^l. 

Let us observe further that while the compactness of F was used in formulation 
of the theorem, it played essentially no role in the proof. This allows us to treat 
other situations, for instance, eigenvalues corresponding to non-straight curves as 
described in Sec. 13.11 Suppose we replace such a F by a family of curves with a 
hiatus, Fg given by the same function 7 where, however, the argument runs over 
K\ (— e, e). Using the above introduced notation, we then have the following result: 

Corollary 3.9. Suppose thatV satisfies the assumptions of Theorem \3.1[ then 
the eigenvalues of Ha,r o-nd Ha^r, obey the asymptotic formula 

Aj(e) = Aj(0) + 2ae sj + o(e) , m{p) < j < n{^) , as e ^ 

Remark 3.10. As before, things look differently in the case of codimension 
two. If we have, for instance, a simple eigenvalue of Ha^r with the eigenfunction 
ip, where F is a curve in MP and perturb the latter by making a 2e-hiatus in it, 
the leading term in the perturbation expansion is again proportional to |<p(0)p, 
however, this time it comes multiplied not by e but rather elne - cf. |EK07| . 

3.5. An isoperimetric problem. The above results do not exhaust ways in 
which the edge shapes influence spectral properties of leaky-graph Hamiltonians. 
Let us mention one more; for simplicity, we restrict ourselves again to the planar 
case, F C K^. If F is of a finite length, we have (Tcss{Ha,r) = [0,oo), while the 
discrete spectrum is nonempty and finite, so that 

(3.6) Ai = Ai(a,F) := inf o-(i/„j) < 0. 

Suppose now that F is a loop of a fixed length and ask which shape of it makes the 
above principal eigenvalue maximal. 

Let us make the assumptions more precise. We suppose that 7 : [0, L] — > M.'^ 
is a closed C^, piecewise smooth curve, 7(0) = j{L); we allow self-intersections 

are other ways to derive such asymptotic expansions, for instance, the technique 
of matching of asymptotic expansions II192| . The advantage of the sketched approach is that it 
requires no self-similarity properties for the family of shrinking sets ■ 
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provided the curve meets itself at a non-zero angle. Furthermore, we introduce 
the equivalence relation: the loops F and F' belong to the same class if one can 
be obtained from the other by a Euclidean transformation of the plane. Spectral 
properties of the corresponding operators Ha.r and Ha^v are obviously the same, 
hence we will speak about a curve F having in mind the corresponding equivalence 
class. The stated assumptions are satisfied, in particular, by the circle, say C := 
{ cos ^ sin ^J^) : s G [0, L] }, and its equivalence class. 

Theorem 3. IF Within the above described class of loops, the principal eigen- 
value Ai(a,F) is for any fixed a > and L > sharply maximized by the circle. 

We will need the following geometric result about means of chords: 

Proposition 3.12. LetT have the properties described above, then 

h{s+u) - 7(s)|f ds < sinP — for pe (0, 2] . 

ttP L 

The right-hand side of the inequality is obviously the value of the integral for the 
circle. Notice that the same is true for loops in R"^ and a similar reverse inequality 
holds for negative powers p G [—2,0). To prove this it is only necessary to establish 
the result for p — 2 which was done in various ways in [Lu66. ACF03., .EHL06] , 
see also |Ex05b| for a local maximum proof. 

Sketched proof of Theorem 13.111 We shall rely again on Proposition l2.3[ 
which relates our eigenvalue problem to the integral equation TZ'^ p0 = on 
L^([0,L]) with TZ'^Y defined similarly as in the proof of Theorem 13.11 we note 
that the operator-valued function k i— > TZ'^ p is strictly decreasing in (0, oo) and 
ll'^^Q rll ^ as K — > cxD. By a positivity improving argument the maximum eigen- 
value of TZ'^ p is simple, and the same is true by Proposition l2 . 31 for the ground state 
of Ha^r- If F is a circle, the latter exhibits rotational symmetry, and using Propo- 
sition 12.31 again we see that the respective eigenfunction of TZ'^^f. corresponding to 
the unit eigenvalue is constant, (^i(s) — L~^^^. Then we have 

maxainllc) ^ i4>i,K]c4>i) ^ 2 ^ K]ci^, s') dsds' , 
while for a general F a simple variational estimate gives 

maxa(7^^V) > i^.K'rM = y / / Tl^r{s,s')dsds' ; 

^ Jo Ja 

hence to check that the circle is a maximizer it sufficient to show that 

Ko{n\T{s)-T{s')\)dsds' > / / KQ{n\C{s)-C{s')\) dsds' 



^0 "'0 "'0 



holds for all K > and F of the considered class. By a simple change of variables 
we find that this is equivalent to positivity of the functional 

L/2 

' Ko{K\r{s + u) -ris)\) - Ko{K\Cis + u) -C{s)\) 



F,(F) :- / du ds 
Jo Jo 
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where the second term is equal to i^of— sin^). Now we employ the (strict) 



convexity of Kq which yields by means of the Jensen inequality the estimate 



^F.(r)> 



rL/2 








'a 



, , , , ^ , , , TTU 
|r(s+u) - r(s)|ds \-Kq \ sin — 



dw , 



where the inequality is sharp unless |r(s+u) — r(s)|ds is independent of s. Finally, 
we note that is decreasing in (0, cxd), hence the result follows from the geometric 
inequality of Proposition 13.121 with p = 1, see |Ex05b| for details. □ 

3.6. Scattering. The investigation of leaky graphs is not exhausted, of course, 
by analysis of their discrete spectrum. Another important problem concerns scat- 
tering on graphs having semi-infinite edges. Our knowledge about this subject is far 
from satisfactory at present and we will concentrate here on a particular situation. 
First of all, we consider again planar graphs, F C M^, only. Secondly, we restrict 
ourself to graphs which can be regarded as a local modification of a straight line. 
And finally, we will analyze the situation which is from the point of view of our 
physical motivation the most interesting, namely the negative part of the spectrum 
where the scattering states are "guided" along the graph edges. 

The unperturbed graph is thus the straight line S = {(a;i,0) : a;i £ K} for 
which s allows a separation of variables; the spectrum is purely a.c. and, in 
particular, the generalized eigenfunctions corresponding to A e are 

(3.7) u;x{xi,X2) = e*(>+"V4)^/^.ie-a|..|/2 

and its complex conjugate Qx. If we perturb the line S locally, we will get a 
nontrivial scattering, but the essentially one-dimensional character of the motion 
will remain asymptotically preserved. Let us specify the perturbation: 

(si ) locality: there is a compact M C such that F \ M = E \ M. 
(s2) finiteness: F \ E is a finite graph with the properties (gl) and (g2). 

We will treat the operator H^x defined by the prescription given in Sec. 12.11 as a 
singular perturbation of -ff^.s supported by the set 

N 

(3.8) A = Ao U Ai with Aq := E \ F , Ai := F \ E = |J F, ; 

i=l 

the coupling constant of the perturbation will take the positive value a on the 
"erased" part Aq and negative one, —a, on added edges Ai. 

For our present purpose we need a more suitable resolvent expression than 
that given by Proposition 12. 3[ instead of (|2.10p we will use the resolvent of i?Q,E 
as the comparison operator. The latter can be expressed, of course, again by 
Proposition 12. 31 we have i?! = Rq -I- ai?^(/ - ai?^^)^^(i?^)*, where for simplicity 
we write i?da;,M = ^'^'^ ^ -"^ "^s is the measure associated with the line by 
A direct calculation then yields the expression 

(T /■ pipx-ip'y I \ 

(3.9) i?|(x-j/) ^ Gfc(x-y) + / ,,,,,, dpAp'^ 



47r3 7r3 (p2_fc2)(p'2__fc2) 2Tfe(pi)- 

with Tkipi) ■— {p\ — fc^)^/^ for any k with Im/c > such that /c^ e C \ [—^0^,00). 
To get the above indicated expression for the resolvent of i/a,r we decompose the 
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Dirac measure associated with the perturbation by p.Sp as 

N 

(3.10) = = + , 

i=l 

where i'q corresponds to the "erased" edge Aq and i^i to Ti . The associated integral 
operator wiU act in the Hilbert space h :— Lp'iv) which inherits from ()3.10p the 
decomposition h = hg © hi with ho := L^{vo) and hi := ®i^iL^{vi). We will 
again need the trace maps, this time of the operator associated with (|3.9p . given by 

(3.11) :h^L2^ = i?| * /i. for /eh 

together with the adjoint (i?|;,^)* : ^ h. and Rj^ ,^1, which is the operator- 
valued matrix in h with the "block elements" R^ij ■ L'^{vj) L'^{vi) defined as 
the appropriated embeddings of p.9p . They have properties analogous to those of 
Proposition 12.31 which can be checked in a similar way, cf. [EK05| . 

Proposition 3.13. The operator Rig ^ is bounded for any k £ (^a, oo). More- 
over, to any a > there is a k^, > such that WRj^^uW holds for n > k„. 

To express the resolvent we introduce an operator- valued matrix in h = ho®hi, 

(3.12) e^--(a~lI + i^|,,J with I:= 

where li are the unit operators in h.^ . By Proposition 13.131 the operator O^'' is 
boundedly invertible for k large enough and we have the following theorem [EK05| : 

Theorem 3.14. Suppose that (6^)^^ e B{\\) hold for fc e C+ and the operator 

R't = R'h + R'hA®^)-\Rl.y 

is defined everywhere in L^(M^). Then belongs to p{Hax) o,nd the resolvent 
{Ha,r — k^)~^ coincides with i?p. 

A simple estimate shows that the operator R^ ^ is Hilbert-Schmidt under our 
assumptions and since the other two factors are bounded by Proposition 13.131 we 
get as a consequence stability of the essential spectrum. 

Corollary 3.15. crcss(i?a.r) = (Toss(i?a,E) = [-ja^,oo) . 

Let us turn now to the proper topic of this section, which is the scattering 
theory for the pair {Ha^r, Ha,T,)- To establish existence of the wave operators by 
the standard Birman-Kuroda method we need to check that the resolvent difference 
B'' := R'h,^{Qiy^{R^,^)* is of the trace class. 

Proposition 3.16. B^'^ is a trace class operator for all k sufficiently large. 

Sketch of the proof. The idea is inspired by |BT92) . We estimate the 
(sign-indefinite) operators (6^")"^ and i?**^ from above and below by a positive 
and negative operator, respectively, which are obtained by taking both signs in 
the matrix multiplying a^^ in (j3.12p the same. Then we can integrate the kernel 
diagonal of the estimates to with a suitable cut-off, using the lemma following 
Thm XI.31 in |RSi79j . and to show subsequently that the trace-class property is 
preserved when the cut-ofi^ is removed, cf. |EK05) for details. □ 
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The existence of wave operators which follows from Proposition 13.161 does not 
tell us much, and we have to find also the on-shell S-matrix relating the incoming 
and outgoing asymptotic solutions. In particular, for scattering in the negative 
part of the spectrum with a fixed A G (—ja^^O) corresponding to the effective 
momentum ka{\) := (A + q;^/4)^/^, the latter are combinations of ujx and ujx given 
by (j3.7p . These generalized eigenf unctions and their analogues ujz for complex 
values of the energy parameter are only locally, of course, but we can use the 
standard trick and approximate them by regularized functions, for instance 

a;f(x) = e~''^icjz(a;) with z e p{Hax) i 

which naturally belong to the domain D{Ha.'si)- Now we are looking for a function 
V'f such that (— Ar — z)^^ = (—As — z)uol. Computing the right-hand side and 
taking the limit limj^o '/'A+ic "^a topology of we find that -0^ still 

belongs to D{Ha,r), and moreover 

where It is the standard embedding from W^'^ to h = L'^{vf^) and R^l'^'^ is the 
integral operator acting on the Hilbert space h, analogous to (|3.1ip . with the kernel 

(3.13) Rl-^^\x-y) lim 4"('+'^)(x-y) ; 

similarly 0'^°('^) := — a^^I — R^^,^ are the operators on h with R^"^^'' being the 
embeddings defined by means of p.l3p . When we remove the regularization, the 
pointwise limit ip\ := lim^^o ^a ceases to be square integrable, however, it still 
belongs locally to and yields the generalized eigenfunction of i?a,r, namely 

(3.14) V-A = C.A + )(0'°*^^)"'^AWA , 

where Jawa is the embedding of lox to L'^{i'\). The on-shell S-niatrix can be 
then found by inspecting the asymptotic behavior of the function ipx as |a;i| — > oo. 
Using the explicit form of the kernel p.l3p derived in |EK02] one arrives by a 
direct computation at the following result: 

Theorem 3.17. For a fixed A G (— ia^,0) the generalized eigenf unctions of 
Ha^r behave under the assumptions (si), (s2) asymptotically as 

where ka{\) '■= (A + a^/4)^/^ is the effective momentum along S and T(A) ,TZ{\) 
are the transmission and reflection amplitudes, respectively, given by 

7^(A) = l-r(A) = ^((eMA))-J..A,JA.A)^. 

4. Strong-coupling asymptotics 

The coupling constant in Ha.r determines how is the particle attracted to 
the graph, and this in turn implies, in particular, which is the "spread" of possible 
eigenf unctions in the direction transverse to the edges. It is thus natural to ask what 
happens in the case of a strong coupling when the wave functions are transversally 
sharply localized. We are going now to show that if the interaction support is a 
sufficiently smooth manifold, various asymptotic formulae can be derived. 
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4.1. Interactions supported by curves. In distinction to the previous con- 
siderations it is rather the dimension than codimension of F that wiU be important. 
As usual, we begin with the case of planar curves, at first finite ones. 

Theorem 4.1. Suppose that 7 : [0, L] ^ is a smooth function, I7I — \, 
which defines a curve V ; then the relation 

ttCTdisc(-ffa,r) = ^ +0(lna) 
zvr 

holds as a — *■ c». In addition, if T is a closed curve without self-intersections, then 
the j-th eigenvalue of the operator Ha behaves asymptotically as 

Xj{a) = — - + Hj + 0{a^^ Ina) , 

where fij is the j-th eigenvalue of the operator Sr ■= ^gp' — jk{s)'^ on L^(0,L) 
with periodic b.c., counted with multiplicity, and fc(s) is the signed curvature ofT. 

Sketch of the proof. Suppose first that F is closed, without self-intersec- 
tions, and consider its strip neighborhood analogous to (|2.7p . in other words, the 
set So onto which the function $0 : [0, L) x (—a, a) defined by 

(4.1) (s,m) 1-^ (7i(s) - w72('S),72(s) +W7j(s)) 

maps, diffeomorphically for all a > small enough. The main idea is to apply to 
Hax the Dirichlet-Neumann bracketing at the boundary of S^, 

(-A^ ) © < Ha.r < (-A^J e Lla, 

where = U A°"' is the exterior domain, and self-adjoint operators 

associated with the forms 



9a*J/] = l|V/||i.(s„)-a J^lfix^dS 

where / G Wq '^(Ea) and W'^''^{T,a) for ±, respectively The exterior \ does 
not contribute to the negative part of the spectrum, so we may consider L^^ only. 

We use the curvilinear coordinates (s,m), the same as in (j2.7p . passing from 
^a.a to unitarily equivalent operators given by quadratic forms 



btjf]^ / / {l + ukis)) 



df 



2 



(s, u) du ds 



ds 

Vis,u)\fis,u){'dsdu-a I |/(s,0)|^ds 





df 


Jo J -a 


du 



2 



(s, u) du ds 



|2, 



with / e W^''^{{0, L) X (—a, a)) satisfying periodic boundary conditions in the 
variable s and Dirichlet b.c. at u = ±a, and 

where V is the usual curvature induced potential |ES89) 



4(l-Hu/c(s))2 2(l-HMfc(s))3 4(H-ufc(s))^ 
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We may employ rougher bounds squeezing Ha,r between H^.a — + T^.a 

with decoupled variables. Here are self-adjoint operators on L?{Q,L) given by 

Ut =-{lTa\\k\U-^^^+V±{s) 

with periodic b.c, where V-{s) < V{s,u) < V^{s) with an 0{a) error, and the 
transverse operators are associated with the forms 

r l/'(")l'dz.-«l/(0)l' 

J —a 

and 

tajf] = tajf] ll^-||oo(|/(a)P + , 

where / € Wg^'^(— a, a) and a, a) for the ± sign, respectively. Their negative 

spectrum can be localized with an exponential precision: there is a c > such that 
has for a large enough a single negative eigenvalue satisfying 

(4.3) - ^ (l + ce— < , < < < (l - 8e--/2) 

To finish the proof, one has to check that the eigenvalues of differ by 0{a) from 
those of the comparison operator, then we choose a = 6a~^ In a as the neighborhood 
width and putting the estimates together we get the eigenvalue asymptotic formula; 
for details see |EY02aj . If F is not closed, the same can be done with the comparison 
operators S^'^ having the appropriate b.c, Dirichlet or Neumann, at the endpoints 
of F; this gives the estimate on 'i (J disc {Ha, r)- D 

The case of a finite curve in M"^ is similar, but we have to be more cautious about 
the regularity of the curve. It will be again a graph of a smooth function, 7 : 
[0, L] M.^ with |7(s)| = 1. To construct the three-dimensional counterpart of the 
"straightening" transformation used in the above proof, we suppose for simplicity 
that F possesses a global Frenet framed and consider the map (pa : [0, L\x Ba ^ 

q^ais, r, e) = 7(s) - r [n(s) cos(0-/3(s)) + b{s) sm{9- f]{s))] , 

where Ba is the circle of radius a centered at the origin and the function /3 has to be 
specified; for small enough a it is a diffeomorphic map on a tubular neighborhood 
Sq of F which does not intersect itself. The geometry of is naturally described 
in terms of the metric tensor gij expressed by means of the curvature k and torsion 
T of F. In particular, in the neighborhood with a circular cross section we can 
always choose the so-called Tang coordinate system, /3 = r, in which the tensor gij 
is diagonal, i.e. the longitudinal and transverse variable decouple |DE95) . 

To state the result we have to note that in the codimension two case the strong 
coupling means large negative values of the parameter a. 

Theorem 4.2. For curves F without self-intersections described above, we have 
tt adisc(i?a,r) = -(-Ca)'/'(l + 0(e"")). 

TT 

as a — > —00, where we put again (^a '■= — 4e^'^^'^"+''^'"'^'-' . //, in addition, T is a 
closed curve, the j-th eigenvalue of the operator Ha,r behaves asymptotically as 

A,(a) = Ca+Mj+0(e""), 



This assumption can be weakened, see the footnote in Sec. 12.41 
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where fij is the j-th eigenvalue of the same operator Sr as in Theorem \4-.l\ 

Sketch of the proof. The argument follows the same scheme. We use 
Dirichlet-Neumann bracketing at the boundary of and estimate the internal 
part using the Tang coordinate system. The effective potential replacing (|4.2p is 
known from [DE95] ; it is important that the torsion does not contribute to its lead- 
ing order as a ^ 0, which is the same as in the two-dimensional case. Also (|4.3p 
has to be replaced by the appropriate two-dimensional estimate, which is again 
exponentially precise, see jEK04aj for details. □ 

The technique used in these proofs can be applied to many other cases. If F is an 
infinite curve, the threshold of the essential spectrum is moved and the estimates on 
tt cTdiscl-ffaT) SLTC uo longer relevant. On the other hand, the eigenvalue asymptotic 
formulaae remain valid under mild additional assumptions [EYOll IEK04a| : 

Theorem 4.3. Suppose that 7 : M ^ R'', = 2,3, satisfies hypotheses of 
Theorems \3.1\ and \3.7[ respectively. In addition, assume that k(s) and k{sY^'^ 
are 0{s~^~^) as \s\ —> 00, and t,t £ L°°(R) for v — Then the asymptotic 
expansions from Theorems and^JT^ hold for all the eigenvalues Xj{a) of H^x! 
when Sr '■= ^gp' ^ jk{s)'^ is now the operator on L^(M) with the domain W^'^(IR). 

Remark 4.4. In this case we need not care about the multiplicity, because 
the spectrum of Sr in L^(R) is simple. This may not be true, of course, in the 
more general case to which the results extend easily, namely for F's consisting of 
disconnected smooth edges, i.e. curves which do not touch or cross each other. 
On the other hand, situation becomes considerably more complicated in presence 
of angles or branchings; we will comment on it in Sec. 17.121 below. 

4.2. Surfaces in R'^. The method works also for interactions supported by 
surfaces, but the geometric part is naturally different. Let us consider first the 
case of a smooth compact and closed Riemann surface F C R"^ of a finite genus 
g. In the usual way |Kli78j . the geometry of F is encoded in the metric tensor 
g^i, and Weingarten tensor /i^ . The eigenvalues k± of the latter are the principal 
curvatures which determine the Gauss curvature K and mean curvature M by 

K ^ detih^'') = k+k_ , M= iTr(/i^'^) = i(fc++fc_). 

The operator H^^r is defined as in Sec. 12.11 For a compact F the essential spectrum 
is [0,cxd) and we ask about the asymptotic behavior of the negative eigenvalues as 
a — > 00. It will be expressed again in terms of a comparison operator: the Sr of 
Theorem 14.11 has to be now replaced by 

(4.4) Sr ■■= -Ar + K - Af^ 

on L^(F, dT), where Ar = —g^^^^d^g^^^g'^'^di, is the Laplace-Beltrami operator on 
F. The j-th eigenvalue /ij of 5'r is bounded from above by that of Ar because 

K = -^{k+ - k^f < 

in analogy with the curve case; in distinction to the latter the two coincide when F 
is a sphere. With these prerequisites we can make the following claim: 
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Theorem 4.5. Under the stated assumptions, jj crdisc(^^Q,r) > j for any fixed 
integer j if a is large enough. The j-th eigenvalue Xj{a) of Ha.T has the expansion 

Xj{a) = ^^'^^ + l^j + 0{a^^ In a) 

as a oo, where fij is the j-th eigenvalue of Sr- Moreover, the counting function 
a I— > tt crdisc(-ffQ.r) behaves asymptotically as 

Ibn 

where \T\ is the Riemann area of the surface T. 

Sketch of the proof. To employ the bracketing, we need to construct a 
family of layer neighborhoods of T. Let {n{x) : x G T} be a field of unit vectors 
normal to the manifold; such a field exists globally because F is orientable. We 
define a map £a : T x (—a, a) — s- R'^ by La(x, u) = x + un(x). Due to the assumed 
smoothness of F it is a diffeomorphism for all a small enough, mapping onto the 
sought layer neighborhood = {a; g : dist(a;,r) < a}. 

By bracketing we get a two-sided estimate for the negative spectrum of Haj^ 
by means of the layer part operators. The latter can be analyzed by means of the 
curvilinear coordinates following [DEKOl] . see |Ex03) for details. One arrives at 
estimates through operators with decoupled variables, ® I -\- 1 ® with 

5*^ := -C±(a)Ar + C^^{a){K - M^) ± va 

and the transverse part which is the same as in the proof of Theorem 14.11 Here 
C±(a) := {1 ± ag~^)'^ with g :— max({||/c+||j^ , ||fc_||^})~-^ and w is a suitable 
constant. The rest of the argument is again analogous to Theorem 14.11 to get the 
counting function one has to employ the appropriate Weyl formula |Ch99] . □ 

Remark 4.6. The connectedness assumption is made for simplicity; the claim 
remains valid if F is a finite disjoint union of smooth compact Riemann surfaces 
of finite genera. Moreover, the asymptotic formula for \\(7d{Hax) is preserved if F 
has a nonempty and smooth boundary, see [Ex03l for a more detailed discussion. 

Under additional assumptions the technique can be applied also to interactions 
supported by infinite surfaces. One possible set of hypotheses looks as follows: 

(asl) injectivity: the map : F x (—a, a) — > 51q C K'^ defined above is 

injective for all a small enough. 
(as2) uniform ellipticity: c^S^^, < g^i, < c+J^i^ for some c± > 0. 
(as3) asymptotic planarity: K, M as the geodesic radius r — ^ oo. 

One can also replace the last requirement by a stronger assumption which implies, 
however, at the same time the validity of (asl): 

(as3') asymptotic direction: the normal vector n ^ uq as the geodesic 
radius r — > oo , where uq is a fixed vector. 

Then one can prove in a similar way as above the following result [EK03j : 

Theorem 4.7. (i) Assume (asl) and (as3), or alternatively (asS'), then we 
have inf (Tess(-ffQ,r) — where e(a) + ja^ = 0(c?%~°"^l'^) as a —> oo. 

(ii) In addition, assume (as2). Unless V is a plane, there is at least one isolated 
eigenvalue of Ha^r below the threshold of the essential spectrum for all a large 
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enough, and moreover, the eigenvalues Xj (a) of Ha,r have the following asymptotic 
expansion, 

Xj{a) = -^a^ + + 0{a^^ \na) 

as a oo, where fXj the j-th eigenvalue of the corresponding operator counted 
with multiplicity. 

In addition to the eigenvalue expansion, we have estabhshed also the existence 
of curvature-induced bound states for non-planar F with suitable spatial asymptotic 
properties in the situation, when the particle is attracted to F sufficiently strongly. 

4.3. Periodic curves. Other large class of leaky graph Hamiltonians for 
which we can investigate the strong-coupling asymptotics in the described way 
concerns periodic manifolds. Let us start with a planar curve F being the graph of 
a smooth function 7 : R ^ M^. In view of the smoothness, the signed curvature 
fc is a function; we assume 

(asl) curvature periodicity: there is L > such that k{s + L) = k{s). 
(as2) curve periodicity: k{s)ds — 0. We may suppose that the nor- 
mal at s is (1,0), then r(- + L) ~ F(-) — (lijh) where period-shift 

components are Ij := sin (^j{2 ~ j) — J*k{u)duj dt. Again without 

loss of generality, we may suppose that li > 0. 
(as3) period cell match: the map (|4.ip is injective for all a small enough 
and $a((0,i) X (-a, a)) C A := (0,/i) x M. 

As usual in a periodic situation we have to perform the Floquct-Blocli decomposi- 
tion [Ku93| . The operator Ha^r{0) on i^(A) is for a 6* e [— 7r,7r) defined through 
the quadratic form as in Sec. 12. 1[ its domain consists of functions u € W^'^{A) sat- 
isfying the boundary conditions u{li, I2 + ■) = e'^u(0, •). In a standard way [EYOl] 
one proves existence of a unitary U : L^(M'^) — > J^_^ L'^i-^) d^* such that 

UH^,tU-^^ / H^,T{S)de and (7(i/„,r) = |J a(i/«,r(0)); 

— 7r,7r) r ^ 

L ' ' [— 7r,7r) 

since 7((0, L)) is compact we have cross(^^Q.r(^)) = [0, 00). Next we need a compar- 
ison operator on the curve analogous to St of Theorem 14. II For a fixed 9 e [— tt, it) 
it is given by the same symbol, 

ST{e) = -^^\k{sf on L^{{Q,L)) 

with the domain {u G W'^'^{{Q,L)) : u{L) = e'''u(0), u'{L) = e'''u'(0)}. We denote 
by ^j{9) the j-th eigenvalue of Sy-{6) counted with multiplicity. Modifying the 
method of Sec. l4.1l to the present situation we get following result [EYOlj 

Theorem 4.8. To any j G N there is an > such that jj crdisc(-ffa,r(^)) > j 
holds for a > aj and any 9 G tt, tt). The j-th eigenvalue of Ha^r{9) counted with 
multiplicity has the asymptotic expansion 

Xj{a,9) ^ -^a^ + fij{9) + 0{a-'^ \na) 

as a 00, where the error term is uniform with respect to 6 ^ [— 7r,7r). 
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Combining this result with Borg's theorem on the inverse problem for Hill's 
equation, we can make a claim about gaps of a{Ha.T)'- 

Corollary 4.9. Assume that T is not a straight line, k ^ Q, then the spectrum 
of Ha,Y contains open gaps for all a large enough. 

In an exactly similar manner one can treat periodic curves in R'^ applying 
the technique to the fiber operator in the Floquet-Bloch decomposition of H^x 
|EK04a| : more care is needed only when defining HaxiP)^ since quadratic forms 
cannot be now used. 

Theorem 4.10. Let T be a periodic curve, without self-intersections and with 
the global Frenet frame, given by a smooth function 7 : R ^ . Suppose, in 
addition, that the period cells Tp of T and A referring to the corresponding operator 
Ha,T match in the sense that Fp = F n A. Then (Tdisc(^^Q,r(^)) has the properties 
analogous to those of the previous theorem, in particular, the j-th eigenvalue of 
Ha,ri^) has the asymptotic expansion of the form 

Xj{a,0) = Ca+ fij{e) +0{e'^°') as a ^ ~oo , 

where fij{9) is the j-th eigenvalue of Sr{d) and the error is uniform w.r.t. 9. 

Remarks 4.11. (i) Corollarv l4.9l has. of course, its three-dimensional analogue. 
The number of gaps which can be open depends on the shape of F. Notice that the 
operator Sr{0)(iO has generically infinitely many open gaps, and in such a 

case the corresponding Ha Y can have any finite number for a large enough. 

(ii) The assumption about a match between the periodic decompositions of the 
curve F and the corresponding may seem restrictive. One should realize, 
however, the period cell need not be rectangular as one usually supposes when 
periodic Schrodinger operators are considered. What we actually need is a complete 
"tiling" of R" by domains with piecewise smooth boundaries. In the case v = 2> such 
"bricks" need not even be simply connected: remember what your grandmother 
was doing with her crotchet to get an example of a curve which is topologically 
inequivalent to a line, or in other words, you cannot disentangle it by any local 
deformation (you can only unwind it by "pulling the ends" ) . 

(iii) The case we have discussed above, namely a single infinite curve in W periodic 
in a given direction, is the simplest possibility. In a similar way one can treat infinite 
families of curves periodic in r < v directions, the only restriction is that their 
components have to satisfy individually the listed assumptions and the distances 
between them must have a uniform positive lower bound. The case r = v differs 
from r < V because then the basic cell is precompact, and therefore the spectrum 
of each i?Q,r(6') is purely discrete. 

(iv) A particular situation occurs when a periodic F consists of disjoint compact 
components. The asymptotic expansions are valid again but now the eigenvalues 
of the fiber comparison operator are independent of the parameter 6 so we have 

and the respective expansion in the three-dimensional case (when the topology may 
be again nontrivial - imagine a chain of interlocked rings). Naturally, the chances 
to have open gaps in this situation are generally better than in the connected case. 

(v) The last comment concerns interpretation of these results. Recall that the 
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deviation of a{Ha.T)i in the negative part at least, from the one described by the 
comparison operator is due to quantum tunnehng. Hence it must be sensitive to 
the appropriate parameter, i.e. the Planck's constant if we reintroduce it into the 
picture. However, the operator — /i^A — v5{x — T) is the multiple of (|l.ip if 
we denote a := vh~'^; in this sense therefore the obtained asymptotic formula 
represent a semiclassical approximation. 

Let us mention one more consequence of these considerations |BDE03] . 

Theorem 4.12. Suppose that the curve T C satisfies the assumptions of 
Theorems \4-8\ and \4-l(^ for v — 2,3, respectively. If v — 1 to any A > there 
is an a\ > Q such the spectrum of the operator Ha.r is absolutely continuous in 
(— oo, —jct^ + A] as long as a > ax. The same is true for v = 3 with —jOi^ replaced 
by Ca provided —a > ax. 

Sketch of the proof. It is easy to check that {Ha,r{(^) ■ G [-tt, tt)} is a 
type A analytic family. The spectral interval in question contains a finite number of 
eigenvalue branches, each is a real analytic function which can be written through 
one of the above asymptotic expansions. The functions fij{-) are nonconstant, hence 
the same is true for Xj{a, •) provided {—l)'^a is large enough. □ 

To appreciate this result recall that if the orbit space of the operator Ha.r 
is compact - cf. Remark I4.11f iv') abov e - there is a wa y to establish the (global) 
absolute continuity of such operators BSSOOi ISSOl) . while the situation with 



r < v is more difficult; we will mention related results in Sec. 16. II below. 

4.4. Periodic surfaces. The treatment of operators Ha.r corresponding to 
periodic surfaces is similar and we describe it only briefly referring to Ex03] for 
details. We consider discrete translations of M.^ generated by an r-tuple {li}, where 
r = 1,2,3. We decompose F, supposed to be a smooth Riemann surface, not 
necessarily connected, and M.^ into period cells Fp and A assuming again that they 
match mutually. The Floquet-Bloch decomposition proceeds as above: we define 
the fiber operators Ha,r{0) on i^(A) through quadratic forms defined on functions 
satisfying the appropriate boundary condition; after that we prove existence of a 
unitary U : L^{R^) /®^ ^y. L'^{A) dO such that 

UHa.rU^'^ = r Ha,riO)de and cr{Ha,r) = [j a{Ha,r{0)) . 

The spectrum of Ha,r{0) is purely discrete if r = 3 while cr ess {H a, r{0)) ~ [0, oo) if 
r ~ 1,2; the eigenvalues are continuous functions of the quasi- momenta 0^. 

As before we need a comparison operator. Its fibers act on L^(Fp,(iF) being 
defined, for instance, by the following prescription 

SriO) g-'/''i-id^ + e^)g^'^ g^'' {-id, + 9,) + K - 

with the domain consisting of e T4^^'^(Fp) such that Ar4> e i^(Fp, dT) satisfying 
periodic b.c. Since Fp is precompact and the curvatures involved are bounded, 
the spectrum of Sr{d) is purely discrete for each 9 G [—tTjTtY; we denote the j-th 
eigenvalue, counting multiplicity, as fJ.j{9). In the same way as above we get 

Theorem 4.13. Under the stated assumptions the following claims are valid: 
(a) Fix A as an arbitrary number if r — 3 and a non-positive one for r — 1,2. To 
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any j G N there is aj > such that Ha,r{0) has at least j eigenvalues below A for 
any a > aj and 9 G [— 7r,7r)''. The j-th eigenvalue Xj(a,9) has then the expansion 

(4.5) Aj(a,6i) ^ ^^a^ + fij{9) + 0{a-Hna) 
as a —> oo, where the error term is uniform with respect to 9. 

(b) If the set cr(S') := [Jg^^_^ a{Sr{9)) has a gap separating a pair of bands, 
then the same is true for a[Ha;r) provided a is large enough. 

4.5. Magnetic loops. Up to now we considered systems witliout external 
fields, which leaves out many situation of experimental interest. For instance, one 
of the often studied features of mesoscopic systems are persistent currents in rings 
threaded by a magnetic flux. For a charged particle (an electron) confined to a loop 
F the effect is manifested by the dependence of the corresponding eigenvalues A„ 
on the flux (p through the loop, conventionally measured in the units of flux quanta, 
2Trhc\e\~^. The derivative dXn/d<j> equals ~\lm where /„ is the persistent current 
in the n-th state. For the ideal loop, e.g., the eigenvalues in absence of other than 
magnetic potential are proportional to (n + 0)^ so the currents depend linearly on 
the applied field. The question is what can we say when the confinement is of the 
type discussed here realized through an attractive interaction on the loop. 

We add a homogeneous magnetic field with vector potential A = ii3(— a;2,xi) 
to our considerations and investigate the Hamiltonian formally given by 

(4.6) Ho,,v{B):^{-iV - Af -a5{x-T) 

in L^(M^). To define it properly we use quadratic form analogous to 



-id^ + ]^By ) ^ 



-idy - -Bx ] ip 



\{I„,i^){x)\^ dx 



with the domain H^^'^(R^). It is straightforward to check that the form is closed and 
below bounded; we identify the self-adjoint operator associated to it with Ha.riB). 

We use the same technique based on bracketing and estimating the operator in 
the strip in suitable coordinates. We need again a comparison operator, this time 

SAB) = - A _ hisf 

on L^{0,L) with V(L-) = e'-^l'^IV'(0+), ^'(i-) = e'-^l"' V'(0+), where Q is the 
area encircled by F. Using it we can state the following result which establishes, in 
particular, the existence of persistent currents on a leaky loop for a large enough. 

Theorem 4.14. Let T be a C^-smooth curve without self-intersections. For a 
fixed j G N and a compact interval I we have ^a{Ha.T{B)) > j for B ^ I if a is 
large enough, and the j-th eigenvalue behaves as 

Xj{a,B) = -^a^ + Hj{B) + 0{a-^ Ina) , 

where Hj{B) is the j-th eigenvalue of Sr{B) and the error term is uniform in B. In 
particular, for a fixed j and a large enough the function Xj (a, •) cannot be constant. 

Sketch of the proof. The argument is closely similar to the analysis of 
fiber operators in Theorem 14.81 the magnetic flux replacing Floquet parameter, 
with small technical differences for which we refer to |EY02bj . □ 
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5. A discrete analogue 

It is often useful to investigate in parallel with leaky graphs analogous discrete 
structures in which the attractive interaction is supported by suitably arranged 
families of point interactions. Let us briefly recall the basic notions, for more 
information and a rich bibliography we refer to (AGHH04) . Consider a selI3 of 
points Y = {yn}nei C M.'^, = 2, 3; if / is infinite we suppose that Y can accumulate 
only at infinity. The operators of interest are point-interaction Hamiltonians Ha,Y, 
typically with the same interaction "strength" at each point, which are defined by 
means of the boundary conditions 

(5.1) Li{i(;,yj) - aLo{ip,yj) ^0, jel, 
expressed in terms of the generalized boundary values 

-^o(V',y):=, lim — -,Li(V',y):= lim [ijj{x) - Lo{ip,y) (j)d{x - y)] , 

where 4>d are the appropriate fundamental solutions, namely 4'2{x) = ln|x| 
and 03 (x) = (47r|a;|)~^ related to the free Green's functions (|2.10p and (|2.13p . 
respectively. The resolvent of Ha^y is given by Krein's formula, 

(5.2) (-i7„,y - fc2)-i = Gk+ J2 [ra,y(fc)]7,' {Gk{-yr), •) Gk{-~yj) 

for fc^ e p{Ha^Y) with Imfc > 0, where Ta.yik) is a closed operator (which is 
bounded in our case) on the matrix representation of which is 

ra,y(fc) \{a - £,1)5,,' ~ Gu{y, - %'(1 - , , , 

where is the regularized Greens's function, ^| — (in ^ — '0(1)) and = 
If a is independent of j, the map k i— > Ta.vik) is analytic in the open upper 
halfplane. Moreover, Ta^vik) is boundedly invertible for Imfc > large enough, 
while for k g C"*" not too far from the real axis it may have a nontrivial null-space. 
By (|5.2p the latter determines the spectrum of the original operator Ha,Y on the 
negative halfline in view the following result analogous to Proposition [2T3] 

Proposition 5.1. (i) A point -n^ < belongs to p{Ha,Y) iff kerPa ^ = {0}. 

(ii) If the operator-valued function k t-^ PQ,,i'(iK)~^ has bounded values in an open 
interval J C with the exception of a point kq G J, where dimkerrQ,_y(iK) = n, 
then — Kq is an isolated eigenvalue of Ha.y of multiplicity n. 

(Hi) an eigenfunction of Ha^Y associated with such an eigenvalue —Kq is equal to 
ij} = '^jGiKoi' ~ yj)y where d = {dj} solves the equation Ta,Y{iKo)d = 0. 

Let us now review discrete analogues of some results derived above. 

5.1. Curved polymers. The operator Ha y referring to a straight equidis- 
tant array Y is called a polymer model in |AGHH04 | . Let us look what happens 
if we abandon the straightness. We adopt the following hypotheses: 

(ad) analogue of (al)-(a3) : Let Y ~ {yn}nez C M", = 2,3, be such 
that \yj — yj+i\ — £ for some £ > 0. It implies \yj — yj' \ < £\j — we 
suppose that the inequality is sharp for some j, j' G Z. Next we assume 



^■^For simplicity we use the symbol Y both for the map / R" and its range. 
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that there is a ci e (0, 1) such that \yj — yj' \ > ci£ \ j ~ and moreover, 
that there are C2 > 0, /i > i, and uj G (0, 1) such that the inequahty 

holds if belongs to the sector Suj of assumption (a2). 

Recall first known facts |AGHH04j about a straight polymer, \yj — yji\ = 
^ I J ~ j'\ for Jt]' ^ ^- Its spectrum is purely absolutely continuous and consists 
of two bands which may overlap if a is not large enough negative. Its threshold 
E^'^ is always negative; in the three-dimensional case it is known explicitly, 

r / , \ -I 2 

while for i/ = 2 we have £2'^ — — ^, where solves the equation 

« + ^(7~ln2) = .g«(0) 

with.g,(e) ^lim^._ {eI-^ i [(" + f (|f)y'-lniv}. 

Theorem 5.2. LetY satisfy the assumptions (ad), then acss{Ha,Y) is the same 
as for the corresponding straight polymer, and the operator Ha,Y has at least one 
isolated eigenvalue below E"'^ for any a E R. 

Sketch of the proof. The argument is a direct counterpart of that used in 
the proof of Theorem l3.ll One checks that the perturbation is sign-definite, pushing 
the spectrum of down, and compact; the result then follows by continuity and 
Proposition cf. [ExOlj for more details. □ 

5.2. Isoperimetric problem. Let us mention also a discrete analogue of the 
problem discussed in Sec. 13.51 we will think now of the point interactions as of 
beads on a loop-shaped string. To be precise, suppose the curve F is the range 
of a function 7 : [0, L] M." which is continuous, piecewise and such that 
7(0) = j{L), and furthermore, |7(s)| = 1 holds for any s £ [0,L] for which 7(5) 
exists. We consider the set Y = {yj} given by 

(5.3) y, ■■^^[^) ' J=0,l,...,iV-l, 

with the indices regarded as integers, yj = ^^(modTV)- A distinguished element of 
the described class is a regular polygon Vn for which the points yj lie in a plane 
C (this is trivial if = 2) at a circle of radius (2 sin ^ . 

We will suppose that the Hamiltonian Ha^y has the property analogous to 
(|3.6p . namely that it possesses a negative principal eigenvalue Xi(a,Y); this is 
automatically satisfied if = 2 while for = 3 it is true if —a is large enough 
}AGHH04] . A counterpart to Theor em 13.111 now reads 

Theorem 5.3. Under the stated assumptions the eigenvalue Ai(a, Yr) is for a 
fixed a and L > globally sharply maximized by a regular polygon, T = Vn- 

We need a geometric result analogous to Proposition 13 . 1 2l which was proved in 
different ways in |Lu66|, IEx06| , see also |Ex05a| for a local proof. 
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Proposition 5.4. Let the set Y be such that \yj+i ~ yj\ < ^, then for any k 
and p e (0, 2] the following inequality is valid 



Notice that the right-hand side is now the value of the sum for the regular 
polygon, and that a similar reverse inequality holds for negative powers p £ [—2, 0). 

Sketched proof of Theorem 15.31 The argument is analogous to that in 
the proof of Theorem 13. Ill All the elements are in place: the ground state is non- 
degenerate and for the regular polygon it has a symmetry, this time with respect to 
a discrete group of rotations which implies that the corresponding eigenfunction of 
ro!,y(«K) is N~^^'^{1, . . . , 1). Using it to make a variational estimate and employing 
the strict convexity and monotony of the resolvent kernel we find that the chord- 
sum inequality of Proposition 15.41 has to be valid with p = 1; a, detailed account of 
the proof can be found in [ExOSaj . □ 

5.3. Approximation by point interactions. Discrete "leaky graphs" de- 
scribed here are useful not only as mathematical objects analogous to the main 
topic of this review. As we are going to mention now, suitable families of them can 
be used to approximate the "true" leaky graphs. Importance of such an approxi- 
mation stems from the fact that apart of particular cases where symmetry allows 
for separation of variables we have no efficient method to find spectral properties of 
the operator Hax- It is true that Proposition 12.31 makes it possible to rephrase the 
original PDE problem as solution of an integral equation of Birman-Schwinger type 
but this is in general a task which not easy either; the purpose of the approximation 
is to convert it into an essentially algebraic problem. 

To get an idea how to proceed in constructing the approximation one can 
compare the spectra of Hax corresponding to a straight line T to that of a straight 
polymer mentioned in Sec. l5.1l above. We let the spacing between point interactions 
go to zero. If the two spectra should coincide in the limit the coupling parameter 
must be inversely proportional to the spacing. It looks queer at a glance but one 
has to keep in mind that the coupling described by the boundary conditions ()5.ip 
becomes weaker as a increases. We have the following result: 

Theorem 5.5. LetT C K.^ be a finite graph obeying the assumptions (gl), (g2) 
and a > 0. Choose k with Imfc > such that the equation a — aR^^a = aR%^ 
has for any ip G L^(R^) a unique solution a which has a bounded and continuous 
representative on T. Suppose next that there is a family {l^jj^i of non-empty 
finite subsets ofT such that \Yn\ :— jj Y,! — > oo and the following relations hold 




N 




for any bounded continuous function f : T 



C, and furthermore 





^ a(y)Gfe(a;-y)-(i?L.„a)(x) ^ 
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as n — > cxD, where \T\ is the sum of all the edge lengths in T. Then the family of the 
operators Ha^.y^ with an '■— |yn|(Q!|r|)~^ approximates the leaky-graph Hamilton- 
ian, Ha^y^ — > H^.v the strong resolvent sense as n — s- cx) . 

Sketch of the proof. The argument is straightforward even if executing it 
needs some effort. We have on one hand the resolvent of Hax given by Propo- 
sition 12.31 on the other hand the resolvent of the approximating point-interaction 
Hamiltonians given by 1)5.2^ : one has to show that their difference applied to any 
ip e L'^{R^) tends to zero as n —^ oo, see [EN03^ for details. □ 

Remarks 5.6. (i) A similar result holds for approximations of leaky surfaces in 
M'^ by families of three-dimensional point interaction, cf. the paper |BFT98] where 
such an approximation was studied for the first time. Both the two-dimensional 
and three-dimensional results are valid for more general sets F and couplings which 
are non-constant functions over such a F. 

(ii) Validity of the theorem extends to two-dimensional systems exposed to a mag- 
netic field perpendicular to the plane. The field need not be homogeneous but it 
has to be sufficiently smooth |Oz06j . 

(iii) As a mean for numerical calculations the above theorem is not optimal. The 
convergence of eigenvalues is slow, roughly as 0{n^^/^), and only in the strong re- 
solvent sense. These flaws can be overcome by considering the family +i?Q„,y„ 
which converges to e^A^ -I- Ha,r in the norm-resolvent sense, and the limit approx- 
imates in the same sense as e ^ [BOOTj. In addition, the numerical 
approximations obtained by this method converge faster. 

5.4. Edge currents in the absence of edges. Magnetic system exhibit 
interesting transport properties manifested, in particular, through edge currents 
discovered in [Ha82, MD S84) and studied in numerous subsequent papers. The 
effect is very robust and can be observed also in situations when the "edge" is thin 
indeed consisting of just an array of point interactions. Consider the situation when 
the interaction sites are arranged equidistantly along a line, for which we can with- 
out loss of generality take the x-axis in 9 a; {xi,X2), and a charged particle 
interacting with them is exposed to a homogeneous magnetic field B perpendicular 
to the plane. In such a situation it is natural to use the Landau gauge in which the 
Hamiltonian can be formally written as 

i^id,, + Bx2f -dl,+Y. -J^) > 

j 

where i is the interaction sites spacing; we write aj to stress that this formal con- 
stant is not identical with the "true" coupling parameter which enters the boundary 
conditions analogous to (jS.ip in a proper definition of the operator. 

The first thing to observe is that the infinitely degenerate eigenvalues which 
the Hamiltonian has in the absence of the point interactions are preserved in the 
spectrum, because one can construct eigenfunctions vanishing at {x^i'^+ji, 0), j € Z. 
Since the system is i periodic, one can perform Floquet decomposition and analyze 
the fiber operators on the strip, which plays here the role of period cell, with a 
single point interaction by Krein's formula analogous to (|5.2|) . arriving thus at the 
following conclusion [EJK99) : 

Theorem 5.7. For any a € R the spectrum of the indicated operator consists 
of the Landau levels B{2n + 1), n = 0, 1, 2, . . . , and absolutely continuous spectral 
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bands; between each adjacent Landau levels there is one such band and the lowest 
one lies below B, the unperturbed spectral threshold. 

Moreover, one can compute the probability current associated with the gen- 
erahzed eigenfunction for a fixed value of the quasi-momentum to see that it is a 
nontrivial vector field describing transport along the array [EJK99 . If the point 
interaction are arranged along a non-straight line, an explicit solution is no longer 
possible but the effect persists - see, e.g., [ChE03] for a regular polygon arrange- 
ment. This is interesting in connection with the result mentioned in Remark l5.6r ii): 
arranging point interactions densely around the loop we can approximate the oper- 
ators discussed in Sec. 14. 5( the "edge" currents are then nothing but the persistent 
currents considered there in the particular case of a strong coupling. 

6. Other results 

6.1. Periodically modulated wires. In Remark I2.ir i) we have mentioned 
that the definition extends to the situation with a non-constant coupling referring 
to the formal expression (jl.ip . Similarly we proceed in the case of codimension 
two; one has to replace a in the boundary condition of (|2.12p by a{s). We will 
denote such operators again Ha.r where we have now a : M — > (0, oo) for 1^ — 2 
and a : M ^ R for = 3. The following result is of a particular intereslF^: 

Theorem 6.1. Suppose that a S L°°(]R) is a periodic function, then the spec- 
trum of Ha r *s purely absolutely continuous; its negative part is non-empty and 
consists of at most finite number of bands. 

Sketch of the proof. The argument is based of investigation of the scat- 
tering for the pair (iJar,— A). Using the Floquet decomposition one establishes 
the existence of wave operators for the fibers and the limiting absorption principle; 
the method employs complexification of the quasi-momentum a la Thomas. To fill 
the details in the case i/ = 2 one has to proceed as in - note 

that the analogous result for a regular potential "ditch" was derived in |FK104| - 
the full proof for = 3 can be found in [EFOTc] . □ 

Let us mention also that in the case v — S one has an analogous result in the 
situation when a periodic interaction is supported by an infinite family of parallel 
lines arranged equidistantly in a plane C K"^ , cf. |EF07c| . 

6.2. A line— and— points model. The number of explicitly solvable models 
in this area is not large, in particular, if we exclude those which can be treated 
by separation of variables. For instance. Theorems 13.11 and 13.171 tell us about the 
discrete spectrum and scattering due to a local perturbation of a straight leaky wire 
but it is difficult to find the eigenvalues or the on-shell S-matrix for a particular 
shape of the deformation. One can achieve more in a caricature model in which 
a straight line is perturbed just by a finite family of point interactions, so the 
Hamiltonian in L^(IR^) can be formally written as 

n 



Recall that for a non-straight periodic F and a constant a we proved the absolute continuity 
only at the bottom of the spectrum provided the interaction is strong enough, cf. Theorem 14. 121 
^^For a recent more general result in arbitrary dimension see | FS07| . 
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where S := {(a;i,0); xi e R}, a > 0, and Y :— {y^'^} is the perturbation support; 
we use again tilde to stress that the fii are not the "true" couphng parameters 
appearing in the boundary conditions analogous to (|5.1|) . 

To analyze such an operator, denoted -ffa,/3, properly defined, we need its re- 
solvent. It can be expressed in a way similar to Proposition 12.31 with the auxiliary 
Hilbert space being now L^(R)©C", or alternatively in analogy with Theorem l3.141 
i.e. by Krein's formula as a rank n perturbation to the resolvent of we refrain 

from stating the explicit formulae which can be found with proofs in [EK04b . Due 
to the finite rank of the perturbation the analysis reduces to an essentially algebraic 
problem. This allows us to prove various results, in particular 

Theorem 6.2. Let f3 = (/3i, /3„) e M" and a > 0. The operator Ha,p has 
a non-empty discrete spectrum, 1 < '^fj{Hajj) < n; the number of eigenvalues is 
exactly n if all the numbers —(3i are large enough. In particular, if n — \ and 
y — (0, a), there is a single eigenvalue; it tends to ~jct^ as \a\ — > oo when := 
_4g2(-27r/3+v(i)) g (_1q,2^o) and to CfJ in the opposite case. 

Recall the usual definition of a resonance as a pole of the analytical continuation 
of the resolvent across (a part of) the continuous spectrum to the lower half-plane. 

Theorem 6.3. Let again n — 1 and y — (0, a), and suppose that > —jo^. 
Then for all \a\ large enough the Hamiltonian has a unique resonance, z{a) = 
fj,{a) — iv^a) with f^a) > 0, which in the limit \a\ oo behaves as 

We will not give details of the proofs referring to [EK04bj for a full exposi- 
tion. In this paper also other spectral and scattering properties of this system are 
discussed as well as an extension to the three-dimensional situation with the inter- 
action supported by a plane and a finite number of points. Furthermore, resonances 
in the case n > 1 are discussed in |EIK07 where one analyzes also the decay of an 
eigenstate of -ffo,/3 due to the presence of the leaky line - see also Sec. 17.191 

6.3. Numerical results. The emphasis in this review is on analytic ap- 
proaches to the problem, hence numerical methods will be mentioned only briefiy. 
As we have mentioned the cases in which the spectral problem can be solved by 
separation of variables are rare and in a sense trivial. Sometimes, however, one can 
divide W into regions exterior to T in which such a separation is possible and to 
find the solutions by the method known in physics as "mode matching" . Consider, 
e.g., a circular T with an angle-dependent coupling; one can use inside and outside 
the circle the Ansatze X^mez '-'^''/n^ ''('') g™*^ where fln^ are suitable Bessel func- 
tions and to find the coefficients Cm^ using the boundary conditions which couple 
the two parts of the plane - see |ET04| for this and similar examples. 

The most versatile method, however, is based on the approximation of leaky- 
graph Hamiltonians by point interactions as we have discussed in Sec. 15. 31 examples 
of spectra obtained in this way can be found in [ENOll IEN031 IOz06| . A case of 
particular interest concerns resonances for Ha,T when F is a straight line with a 
buckling having a narrow bottleneck. Of course, the approximation method applies 
to finite F's only but one can use the approach popular among the physicists — 
rigorously justified so far only for one-dimensional potential scattering [HMOO] 
— in which one cuts the system to a finite length and observes the dependence 
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of eigenvalues on the cut-off size; resonances are manifested by a pattern in which 
intervals of almost constancy are interlaced with steep jumps. In our case resonances 
appear when the bottleneck half- width becomes comparable with the characteristic 
transverse size of the generalized eigenfunction (|3.7p . i.e. {2a)~^, cf. |EN03j . 

7. Open problems 

While we have reviewed numerous results in the previous sections, the topic 
we are discussing here is new and many questions remain still open. In this closing 
section we are going to list some problems the reader may wish to learn about, 
or better, to attack. As usual in mathematics, the number of generalizations is 
unlimited and we restrict only to those which we regard as both meaningful and 
reasonably close to the above exposition. Even in such a class one can find very 
different problems. Some of those listed below are rather technical and require 
mostly perspiration to achieve the result; others are more challenging when we do 
not know what the result might be, which method to apply, or possibly neither 
of these - in short, one cannot start without a proper inspiration. Presenting the 
problems, however, we will not rank them according to these criteria but rather list 
them in the order they appear in the text. 

7.1. Approximation by regular potential "ditches". This is an impor- 
tant link between our singular model and description of realistic systems of "wires" 
of a small but finite widtlF^. In Theorem 12.21 we considered the simplest case of 
a single curve F C M'^. If F has angles, or even it is a nontrivial graph with 
branchings, one expects the analogous result to be valid, however, the method used 
in the proof is no longer applicable and it seems likely that the edge and vertex 
parts of the approximating potentials have to be treated separately. Similar results 
can be expected if the "edges" of F are surfaces in R^. 

A related question concerns existence of curvature-induced bound states in 
such "fat" leaky graphs, as well as counterparts of other results discussed in Sec. [3l 
Using the above mentioned convergence result, combined possibly with minimax 
estimates, one can establish the analogue of Theorem 13.11 provided the potential 
"ditches" in question are sufficiently deep and narrow; it is less clear whether the 
claim will remain true generally outside the asymptotic regime. 

On the other hand, a potential approximation to H^.r is more delicate in 
case of codimension two. A hint can be obtained from approximations of a two- 
dimensional point interaction |AGHH04| where one starts from a potential well 
having a zero-energy resonance and scales it in a particular nonlinear way. One can 
conjecture that this would yield an approximation for a single smooth F if we take 
such a potential family in the normal plane to F, while the problem may be more 
complicated in presence of angles and branchings. 

7.2. More singular leaky graphs. The previous problem brings to mind a 
related question. Transverse to the graph edges we choose the S interaction as the 
mean to describe the way in which the particle is attracted to F. This is not the 
only possibility, though, even if we restrict ourselves to interactions supported by a 

^^In the photonic-crystal setting, this problem is discussed in IFKu98| . The relation between 
the leaky-graph Hamiltonian and the corresponding pseudo-differential (Dirichlet-to-Neumann) 
operator is another foundational issue which deserves attention, especially in view of the motivat- 
ing numerical results obtained for particular geometries (KK02] . 
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single point. We know that such singular interactions form a four-parameter family 
containing some prominent cases such as the 5' interaction |AGHH04| . Let us 
recall that the latter are not mere mathematical artefacts as it can be seen, e.g., 
from the fact that they are approximated by triples of properly scaled 5 interactions 
|CS98j and even by regular potentials [ANOOllENZOT] . 

It is appropriate at this place to mention that we have not discussed here 
counterparts of leaky-graph Hamiltonians with a repulsive interaction, formally 
—A + a5{x — r) with a > 0, which can be introduced in the case codimF = 1 
as in Sec. 12.11 the reason is that the model has different spectral and scattering 
properties as well as physical interpretation. A combination of the attractive and 
repulsive interactions on a "triplicated" graph where each edge of F is replaced by 
three close edges approximating in the cut the 5' interaction according to |CS98| 
makes perfect sense, however, and justifies interest in such models. 

For 5' and more general leaky graphs we do not have at our disposal a "natural" 
quadratic form definition analogous to (|2.3p or a Birman-Schwinger-type expression 
for the resolvent, hence it is a priori not clear which ones of the spectral and 
scattering properties discussed here can be extended to this case. 

7.3. More on curvature induced spectra. The existence of a discrete spec- 
trum for an infinite curve F which is not straight but it is asymptotically straight 
in a suitable sense as in Theorem 13. II is an interesting result, however, a full under- 
standing requires more. We already know that jj crdisc(^Q,r) can be made arbitrarily 
large finite for a curve as simple as -ff2(/3) of Sec. 13.21 It is intuitively clear that 
a rich discrete spectrum is to be expected when the edges of F come close to each 
other at numerous places or over long stretches; it would be desirable to have a 
more quantitative expression of this intuitive statement. 

It is useful to stress at this place that the leaky curves have a lot in common 
with quantum waveguides (ES89, DE95J. In both cases we study solutions to the 
Schrodinger equation localized in a tubelike region, here in a "soft" way through 
an attraction to a curve in contrast to the "hard" way with Dirichlet boundary 
conditions. This analogy makes is easier to understand effects like existence of 
localized solutions due to curvature because they are similar in both caseJ^. This 
analogy suggests various questions, for instance, about small bending asymptotic 
behavior. Take a curve which is straight outside a compact and differs only slightly 
from the straight line; an archetype of such a behavior is a broken line corresponding 
to H2{'K — 9). Using the method from the proof of Theorem 13. II one can check that 
for 9 small enough such a Hax has a single eigenvalue \{9). We conjecture that 

(7.1) X{9) = -^a' ~ cO^ + 0{e'') 

in analogy with the slightly curved waveguide [ABGM9TI IDE95| . The positive 
constant c here depends on a, in particular, for H2{'k — 9) we have c = da^. The 
methods discussed in this paper, however, do not give a way to prove the conjectured 
asymptotic behavior and another approach is needed. 



This is not to say that all properties are the same. As example is provided by nodal lines of 
the discrete— spectrum eigenfunctions: in thin waveguides they can be described by means of the 
one-dimensional comparison operator IFKrOT] while for leaky curves they extend over the plane 
and it is more complicated to analyze their behavior. 
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The assumption (a2) in Sec. I3.1l was used to ensure that the curvature induced 
spectrum is discrete. The proof of Theorem [3ll] shows, however, that any departure 
from a straight hne pushes the spectrum dowrFl. the question is which character 
does it have below —ja'^. The cases to be analyzed include, e.g., curves with a slow 
curvature decay or a line with sparse deformations. Analogies with one-dimensional 
Schrodinger operators suggest that one might get different spectral types, however, 
they have to be used with caution since the other dimension may play a rolj^. 

The above question belongs to those which are open also in the quantum waveg- 
uide setting mentioned above. Another problem of this sort concerns the spectral 
multiplicity, specifically, one would like to know whether - or possibly under which 
conditions - is the discrete spectrum of Theorem 13.11 simple. 

7.4. More on star graphs. The example discussed in Sec. 13.21 is simple, 
nevertheless, some open questions remain. Notice that the important parameters 
are the angles only, because a star-shaped T is self-similar and a change of a is 
equivalent to a modification of the length scale. One can ask, for instance, about 
the configuration which maximized the ground state; a natural conjecture is that 
this happens in the case of the maximum symmetry, f3i = 2tt/N, i — 1, . . . , N . Less 
obvious is the answer to the question about existence of closed nodal lines which are 
expected if N is large enough. A numerical example [EN03| based on Theorem 15. 51 
suggests this happens, e.g., with the fourth eigenfunction of iJio(7r/5). One can 
ask what is the minimum N for which a leaky-star Hamiltonian has a eigenfunction 
with a closed nodal line, about stability of nodal patterns w.r.t. the angles, etc. 

7.5. Bound states for curved surfaces. The analogy with quantum wave- 
guides extends to the case of a surface in R'^ : by Theorem 14.71 we know that non- 
planar strongly attractive surfaces satisfying (asl)-(as3), or (as2) and (as3') give 
rise to a non-empty discrete spectrum. Recall that for Dirichlet quantum layers we 
also do not have an universal result; they bind if their width is small enough, if 
the total Gauss curvature is non-positive, or under various symmetry assumptions 
[DEKOll ICEKCM] . The thinness assumption is the analogue of a strong coupling, 
a — > oo, discussed here. One is interested whether and under which conditions does 
the curvature imply existence of bound states for a fixed coupling parameter a. 

Guided again by the quantum layer analogy [EKrOl] one can ask about the 
weak-coupling asymptotics of bound states corresponding to mildly curved sur- 
faces. In distinction to the conjecture (j7.ip one guesses that in the two-dimensional 
case the binding will be exponentially weak with respect to a suitable deformation 
parameter. The method to prove this result is not clear, though. 

7.6. Bound states of nontrivial graphs of codimension two. We have 
seen repeatedly that leaky graphs of codimension two require more subtle analysis. 
This applies already to the definition of i/^.r which we have briefly described in 
Sec. 12.41 for the case when F is a curve. There is little doubt that the boundary 
conditions (|2.12p can be used for a non-trivial graphs F as well, however, the proof 
using the resolvent formula (|2.14|) has to be worked out properly. 



Or alternatively, for a pure mathematician, to the left — the same for the previous text. 
^^Interesting effects can be seen also in the positive part of the spectrum. Recall the situation 
where T is an infinite family of concentric, equidistantly spaced circles, then the cress{Hc,^r) consists 
of interlaced intervals of absolutely continuous and dense pure point spectrum [EF07b] . 



LEAKY QUANTUM GRAPHS 



35 



Another question concerns the possible analogue of Corollary 13.21 Due to the 
absence of a "natural" quadratic form associated with Ha^r it is not obvious that 
bound states are preserved when edges are added to a non-straight curve creating 
a non-trivial graph, although one expects that this will be the case. 

7.7. Geometric perturbations for general graphs. The results discussed 
in Sec. 13.41 are formulated for a single smooth curve or manifold. With an additional 
effort one can extend validity of the asymptotic formula to a general graph T 
provided the point at which the hiatus is centered is in the interior of an edge. 
While a similar claim may be still true for edges "narrowly disconnected" in a 
vertex, a more subtle analysis is needed to see whether such a claim is valid. 

7.8. More on isoperimetric problems. Theorem 13 . 11 1 raises naturally vari- 
ous questions about possible extensions. A modification for a loop in M.^ is relatively 
straightforward and one can ask also whether there are local maxima for loops of 
a fixed knot topology. On the other hand, the corresponding problem for a closed 
surface in M'^ seems to be more complicated and direct generalizations do not work. 
Returning to the two-dimensional situation, and referring to Sec. 17.41 above, one 
can ask more generally about configurations which maximize the ground state for 
r with a fixed topology and edge lengths, etc. 

7.9. More on absolute continuity. The case of with a periodic curve 
r remains to be one of the important challenges in this area. One certainly expects 
that the spectrum of such an operator will be purely absolutely continuous, however, 
proof of this assertion is missing. So far we have only the result in the related but 
different case expressed by Theorem 16. 11 and a partial answer to the question given 
by Theorem l4.121 let us note that the absolute continuity is expected be valid even 
for r consisting of infinitely many disconnected finite components for which the 
strong-coupling argument does not work. 

7.10. More on scattering. The setting we choose in Sec. [321 has the advan- 
tage of simplicity offered by the straight- line comparison operator Ha,^,. Scatter- 
ing for a general leaky graph with straight leads outside a compact can be treated 
similarly, but needs a more elaborate formulation: one has to compare motion 
on each external lead with that of Ha.Yi using a suitable identification operator. 
The situation becomes even more complicated when the leads are not straight but 
only asymptotically straight, in particular, if the "curvature decay" is slow so the 
curvature-induced interaction is of a long range. 

A different class of problems concerns scattering at positive energies. Here one 
has to distinguish two cases. If T is compact we deal with standard scattering 
on a potential, albeit a singular one; this situation was discussed in ^BT92]. On 
the other hand, if there are infinite edges the situation is more complicated due to 
presence of guided states; note that they can exist also at positive energies, just 
take (|3.7p with A > 0. This case remains so far largely untreated. 

7.11. Strong coupling behavior of scattering. Speaking about scattering 
one may ask whether there are asymptotic formulae analogous to those we derived 
for the discrete spectrum in SecIH We again suppose that F is a connected infinite 

smooth curve straight outside a compact. For large a the wave function is 
localized in a small neighborhood of F and one can expect that the reflection and 
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transmission amplitudes which determine the scattering in the negative part of the 
spectrum will be determined in the leading order by the local geometry of the curve. 

Let iSr.a(A) be the on-shell scattering matrix at energy A for the operator Ha^r, 
compared to the free dynamics in the leads, and denote by Ss^ (A) the corresponding 
quantity for the one-dimensional comparison operator of Theorem 14.31 One 
conjectures that for a fixed k ^ and a — s- cx) we have the relation 



to prove it and to find the convergence rate one has to show that the correspond- 
ing generalized eigenfunction of Ha.r and that of Sr multiplied by the transverse 
eigenfunction analogous to that in (|3.7p converge to each other, and how fast. 

7.12. More on strong coupling asymptotics. This brings us to strong- 
coupling asymptotic results of Sec. [H One question left open there concerns man- 
ifolds with boundaries; for simplicity let us consider a finite curve F which is not 
closed. As we have pointed out in the proof of Theorem 14.11 the direct use of the 
bracketing technique gives the first assertion of the theorem but they are not precise 
enough to yield an asymptotic formula. We conjecture that it is again of the form 



where /ij is now the j-th eigenvalue of the operator given by the same expression, 
but with Dirichlet boundary conditions. A way to prove this result would be to 
consider a bracketing on prolonged tubular neighborhoods extending to the distance 
of order a over the endpoints of F; a trouble to overcome is that one cannot then 
separate variables in the leading order. 

All the considerations of Sec. [4] required a sufficient smoothness of F. If this is 
not true the results are no longer valid. As an example consider again the operator 
H2if3) of Sec. 13.21 using a scaling transformation we find that fij in the above 
formula has to be replaced by — (A^ + i) where Xj is the j-th eigenvalue of 
H^iP) corresponding to a = 1. We conjecture that if an otherwise smooth F has 
one angle equal to tt — /?, the j-th eigenvalue asymptotics of the corresponding 
is again Aj-a^ -I- o(q;^); if there are more angles the situation is similar but the 
numbering of the eigenvalues changes in the appropriate way. 

7.13. Strong-coupling graph limit. So far we have considered strong cou- 
pling for a single curve. If F is a nontrivial graph with branchings, one expects a 
behavior similar to that of a curve with angles: each vertex will contribute to the 
asymptotics below —^a^ by XjO^ -l-o(a^) where Aj's are now the eigenvalues of the 
corresponding leaky star-graph Hamiltonian Hn{(3) with n being the vertex degree. 

This leads us to the question, whether in the limit a ^ oo one can get a mean- 
ingful expression of the standard quantum graph type, i.e. an operator on L^(F). 
The question makes sense, of course, only if we perform a suitable energy renormal- 
ization. The most natural way to do that is to stay in the vicinity of the "transverse 
threshold", i.e. to subtract the diverging factor — jo^- If the above conjecture is 
valid the limit will be generically trivial, i.e. a graph with Dirichlet-decoupled 
edges. There could be nontrivial limits, however, in cases when the involved fam- 
ily of leaky-graph Hamiltonians has a threshold resonance (or a singularity which 
stays in the vicinity of such a resonance in the limiting process). The problem is 
analogous to the squeezing of Dirichlet fat graphs mentioned in the introduction, 




Aj [a) 



— + fij + 0{a ^ In a) , 
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cf. [PoOSj, ICE07|, IGr07| . and it is likely to be no less difficult. One can also 
conjecture that if the reference point is chosen instead as with A > — jo^ one 
gets generically a nontrivial limit with the vertex conditions determined by the 
scattering on the appropriate leaky star graph in analogy with [MV07| . 

7.14. More on resonances. The presence or absence of resonances is one of 
the most important features in scattering. In this survey we touched that subject 
but not very deeply: we were able to establish the existence of resonances in the 
caricature model of Sec. 16.21 and we quoted numerical results which suggest that 
resonances are present in some other situations. The first problem here is whether 
one is able to establish the existence of resonances for a wider class of operators 
Ha.T with an infinite F having some number of (asymptotically) straight "leads" . 

The question has to be put more precisely. The resonances in Theorem 16.31 are 
understood in the common sense as poles of the analytically continued resolvent. 
For the model in question one can check easily |EK04bj that they are at the 
same time singularities of the scattering matrix. A similar equivalence between the 
resolvent and scattering resonances is expected to be valid for other leaky-graph 
Hamiltonians as well, and its verification should accompany the existence proof. 

One more problem concerns the method we mentioned in Sec. 16. 3[ called usu- 
ally the -approach by the physicists, which allows to identify resonances by 
spectral methods. As we have noted the method is rigorously justified for the one- 
dimensional potential scattering only [HMOO| and it is naturally desirable to find 
an appropriate formulation and proof of it in the present situation. 

7.15. More on magnetic leaky graphs. This is another subject not much 
explored so far, and also one where the analogy between leaky graphs and quantum 
waveguides sometimes fails, e.g., a waveguide in a homogeneous field has dominantly 
a continuous spectrum due to edge states "skipping" along the boundary while the 
spectrum of its leaky-graph counterpart is a point one being dominated by the 
Landau levels. A numerical example worked out in [ET04j using mode matching 
indicates that the asymptotic behavior described in Theorem 14. 141 can be destroyed 
if we keep a fixed, albeit large, and make the magnetic field strong. 

In connection to that one can ask, e.g., what will happen with the curvature- 
induced bound states of Theorem 13.11 if the system is exposed to a magnetic field 
a la (|4.6p . in particular, whether they will survive an arbitrarily strong field. In 
a similar vein, it is possible to ask what will happen with the curious magnetic 
transport described in Theorem 1 5. 71 if the point interaction array in question is not 
straight but "bent" and only asymptotically straight. 

A much weaker perturbation is represented by local magnetic fields. Here again 
one can derive an inspiration from the waveguide theory: it is known that spectrum 
of a waveguide with such a field is, contrary to the non-magnetic case, stable with 
respect to small perturbation as a consequence of a Hardy- type inequality [EKo05| . 
One can ask whether a similar result is true for leaky curves slightly different from 
a straight line if a local magnetic field, regular or singular, is present. 

7.16. Perturbations of periodic graphs. In addition to the absolute conti- 
nuity problem mentioned above periodic graphs pose other questions. The spectrum 
is expected to have a gap structure, and in some cases one is able to establish ex- 
istence of open gaps as, e.g., in Corollary 14.91 In such a case one can ask about 
the effect of local perturbation to such a periodic F, in particular, under which 
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conditions they give rise to eigenvalues in gaps, what is the number of the latter 
and their dependence on the perturbation parameters, etc. 

7.17. Absence of embedded eigenvalues. Another question which arises in 
connection with the previous problem, but it it is not restricted to perturbations of 
periodic graphs, concerns embedded eigenvalues. It is well know that in the "usual" 
quantum graphs the unique continuation principle in not valid [KV06| . and as a 
consequence, we encounter frequently situations with eigenvalues embedded into the 
continuum and associated with compactly supported eigenfunctions. Leaky graphs 
are different and one naturally asks whether one can establish generally absence of 
such eigenvalues, and also about possible existence examples if there are any. 

7.18. Random leaky graphs. Quantum graphs with various sorts of ran- 
domness have been recently an object of intense interest - see, e.g., |ASW06l 
IEHS071 IKP071 IGV07| and references therein. In contrast, not much was done 
about the analogous problem for leaky graphs which are certainly of interest be- 
cause they mix elements of a one-dimensional and multidimensional behavior. One 
naturally expects a localization if either the edge shapes or the coupling constants 
on them become random, and there is a numerical evidence supporting these ex- 
pectations [ET04] . The truly important question, however, concerns the existence 
and properties of a mobility edge in such systems. 

The effects of randomness can be also studied through a point interaction coun- 
terparts to the operators iJ^.r- If we take a point-interaction polygon with a mag- 
netic field mentioned in Sec. 15.41 and randomize the coupling constant, one can 
observe numerically how the transport is destroyed. One can ask therefore what 
can be proven about the spectrum of the model from Theorem 15.71 if the coupling 
constants are made random, whether there is a localization and whether a part of 
the absolutely continuous spectrum will survive - recall that for a similar model 
with an array replaced by a rectangular point-interaction lattice a localization was 
proved in the low- lying spectral bands |DMP99| . 

7.19. Time evolution. So far we spoke about the stationary aspect of the 
problem, even when we discussed problems such as scattering. We do not know 
much about the way the wave functions evolve in leaky-graph systems. One can 
ask, for instance, about smoothing properties analogous to those of the usual 
Schrodinger operators. In some cases this is true — take the example of Hax 
corresponding to a straight line in the plane — but a caution is needed since it 
is known that the time evolution in systems with singular interactions can have 
sometimes rather counterintuitive properties [EF07a) . 

A more specific question is associated with resonances in quantum graphs. The 
usual duality between resonances and unstable states motivates us to ask what 
happens if we fix a resonant state, typically an eigenstate of an unperturbed oper- 
ator associated with Ha^r embedded into the continuum, at an initial instant, and 
ask about the way in which it decays. For the caricature model of Sec. 16.21 it can 
be done |EIK07j : one would welcome to learn more about the decay processes in 
physically more interesting cases. 

This problem list is in no case exhaustive. If the reader made it to this point 
— my greetings to such a persistent colleague — I am sure he or she managed to 
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formulate many additional questions on the way, and is ready to address them. It 
remains for me only to express a good luck wish in such an effort. 
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